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Quantized linear systems on integer lattices: 
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In loving memory of my father, Gennadiy Ivanovich Vladimirov (21.04.1933 - 24.12.2014) 

Abstract 

The roundoff errors in computer simulations of continuous dynamical systems, which 
are caused by finiteness of machine arithmetic, can lead to qualitative discrepancies be¬ 
tween phase portraits of the resulting spatially discretized systems and the original systems. 
These effects can be modelled on a multidimensional integer lattice by using a dynamical 
system obtained by composing the transition operator of the original system with a quan¬ 
tizer representing the computer discretization procedure. Such model systems manifest 
pseudorandomness which can be studied using a rigorous probability theoretic approach. 

To this end, the lattice is endowed with a class of frequency measurable subsets and a 
spatial frequency functional as a finitely additive probability measure describing the limit 
fractions of such sets in large rectangular fragments of the lattice. Using a multivariate ver¬ 
sion of WeyTs equidistribution criterion and a related nonresonance condition, we introduce 
an algebra of frequency measurable quasiperiodic subsets of the lattice. The frequency- 
based approach is applied to quantized linear systems with the transition operator Ro L, 
where L is a nonsingular matrix of the original linear system in and R : ^ 17^ 

is a quantizer (in an idealized fixed-point arithmetic with no overflow) which commutes 
with the additive group of translations of the lattice. It is shown that, for almost every 
L, the events associated with the deviation of trajectories of the quantized and original 
systems are frequency measurable quasiperiodic subsets of the lattice whose frequencies 
are amenable to computation involving geometric probabilities on finite-dimensional tori. 
Using the skew products of measure preserving toral automorphisms, we prove mutual in¬ 
dependence and uniform distribution of the quantization errors and investigate statistical 
properties of invertibility loss for the quantized linear system, thus extending the results 
obtained by V.V.Voevodin in the 1960s. In the case when L is similar to an orthogonal 
matrix, we establish a functional central limit theorem for the deviations of trajectories 
of the quantized and original systems. As an example, these results are demonstrated for 
rounded-off planar rotations. 


*This work is a slightly edited version of two research reports: I.Vladimirov, “Quantized linear systems on 
integer lattices: frequency-based approach”. Parts I, II, Centre for Applied Dynamical Systems and Environmental 
Modelling, CADSEM Reports 96-032, 96-033, October 1996, Deakin University, Geelong, Victoria, Australia, 
which were issued while the author was with the Institute for Information Transmission Problems, the Russian 
Academy of Sciences, Moscow, 127994, GSP^, Bolshoi Karetny Lane 19. None of the original results have been 
removed, nor have new results been added in the present version except for a numerical example on p. 58 of the 
last section. 

fPresent address: UNSW Canberra, ACT 2600, Australia. E-mail: igor.g.vladimirov@gmail.com. 
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1 Introduction 

The roundoff errors in computer simulations of continuous dynamical systems, which are caused 
by finiteness of machine arithmetic, can lead to dramatic discrepancies between phase portraits 
of the resulting spatially discrete systems and the original systems ^ IH dH. A model class 
of these spatially discretized systems is provided by dynamical systems on multidimensional 
integer lattices, obtained by composing the transition operator of the original system with a 
quantizer which represents the computer discretization procedure flU [T3l [T71. Trajectories of 
such model systems exhibit pseudorandomness. This suggests that their statistical properties 
can be studied from a probability theoretic point of view by equipping the integer lattice with 
a probability measure so as to quantify the discretization effects for a sufficiently wide class of 
dynamical systems. 

The present study is concerned with those “events” on the lattice which are associated with 
the deviation of trajectories of the discretized and original systems. We develop a rigorous 
probabilistic approach to quantized linear systems on the n-dimensional integer lattice Z” with 
the transition operator R o L, where L is a nonsingular matrix of the original linear system in 
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and the map i? is a quantizer in an idealized fixed-point arithmetie with no overflow. To 
this end, the lattice Z" is endowed with a class of frequency measurable subsets and a spatial 
frequency functional as a finitely additive probability measure on such sets. We introduce an 
algebra of quasiperiodic subsets of the lattice, which are frequency measurable under a non¬ 
resonance condition, and apply this frequency-based approach to the quantized linear systems. 
We show that almost every matrix L satisfies the nonresonance condition, and the events asso¬ 
ciated with the deviation of trajectories of the quantized and original systems over a bounded 
time interval are representable as frequency measurable quasiperiodic subsets of the lattice. The 
frequencies of these sets are amenable to computation which involves geometric probabilities 
on finite-dimensional tori. In particular, for the generic matrices L, satisfying the nonresonance 
condition, we prove the mutual independence and uniform distribution of quantization errors, 
and investigate the statistical properties of the invertibility loss for the transition operator. In the 
case where L is similar to an orthogonal matrix, this allows a functional central limit theorem 
to be established for the deviations of trajectories of the quantized and original systems. As 
an illustration, we apply these results to a two-dimensional quantized linear system associated 
with the problem of rounded-off planar rotations. 

Given the length of this report, we will now outline its structure and main results. In Sec¬ 
tion |2l the n-dimensional integer lattice Z"^ is endowed with di frequency functional F. The 
frequency F(A) of a set A C Z"^ is defined as the limit fraction of points of the set in unbound¬ 
edly increasing rectangular fragments of the lattice. The class S of frequency measurable sets, 
for which this limit exists, is closed under the union of disjoint sets and the complement of a 
set to the lattice, and F is di finitely additive probability measure on the class S os such subsets 
of the lattice. The frequency F is defined in Section 12.21 as an average value of the indicator 
function of a set, with the average value functional A on a class of averageable functions being 
introduced in Section |2T] Since S is not a cr-algebra, and F is not a countably additive measure, 
the triple (Z”, S, F) can be regarded as an unusual probability space (where the average value 
functional A plays the role of expectation) which does not satisfy Kolmogorov’s axioms [fT4]| . 
Nevertheless, this probability space allows a random element with values in a metric space X 
to be defined as a map : Z*^ —)■ X for which there exists a countably additive probability mea¬ 
sure D such that the preimage of every D-continuous (jU Borel set i? C X is frequency 

measurable with frequency = D{B). Any such map g is called D-distributed and 

generates an algebra C iS of events related to the behaviour of g. The distributed maps 
are discussed in Section 12.31 The averagability of a function, the frequency measurability of 
a set and the property of a map to be distributed are closely related to each other and can be 
formulated in terms of the weak convergence of probability measures in metric spaces [[T] . 

Section |3] introduces a class of quasiperiodic objects which are nontrivial examples of av¬ 
erageable functions, frequency measurable sets and distributed maps. Section describes 
quasiperiodic subsets of the lattice TP. Any such set is specified by a Jordan measurable set 
G C [0,1)™ and a matrix A G and is denoted by Qm{G, A). The set Qm{G, A) is formed 

by those points a; G Z” for which the m-dimensional vector Ax, whose entries are considered 
modulo one, belongs to G. The resulting class Qm{A) of A-quasiperiodic sets is an algebra 
with respect to set theoretic operations. Furthermore, Section [33] shows that if the matrix A is 


nonresonant in the sense that the rows of the matrix 


In 

A 


are rationally independent (with the 


identity matrix of order n), then every A-quasiperiodic set Qm{G, A) is frequency measurable 
and its frequency coincides with the m-dimensional Lebesgue measure meSmC of the set G. 
Similarly, a A-quasiperiodic function on Z"^ is defined as the composition / o A of the linear 
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map, specified by a matrix A G and a mesm-continuous bounded function / : —)■ M, 

unit periodic in its m variables. In a similar vein, Section [A4l defines a A-quasiperiodic map 
from to a metric space X as the composition o A of the linear map A with a unit peri¬ 
odic mesrn-continuous map g : M'” —)■ X. It is shown that if the matrix A is nonresonant, 
then every A-quasiperiodic function / o A is averageable, with A(/ o A) = f{u)(lu, 

and any A-quasiperiodic map o A is distributed, with the algebra SgoK consisting of of A- 
quasiperiodic sets. It turns out that mesmn-almost all matrices A G are nonresonant, and 
the corresponding algebras Qm(A) of A-quasiperiodic sets are parameterized by Jordan mea¬ 
surable subsets of [0,1)™. Therefore, a typical matrix A leads to sufficiently rich classes of 
A-quasiperiodic frequency measurable sets, averageable functions and distributed maps, which 
can be studied in the framework of the probability space (Z"^, 5, F). The above properties of 
quasiperiodic objects are established by using a multivariate version of what is known as the 
method of trigonometric sums in number theory [fT^ and Weyl’s equidistribution criterion in 
the theory of weak convergence of probability measures lUl. The quasiperiodic sets are then 
extended to infinite dimensional matrices A. More precisely, for a matrix C G the class 

Qoo{C,) of /^-quasiperiodic sets is defined as the union of the algebras Qrn(A) over all positive 
integers m and all submatrices A G of the matrix C. The class Qoo(/^) is also an al¬ 
gebra and consists of frequency measurable subsets of the lattice Z”, provided C G is 

nonresonant in the sense that all its submatrices A G are nonresonant. The discussion 

of quasiperiodic objects in Sections [T2l43.41 employs the notion of a cell in M™, which extends 
that of a space-filling poly tope [I3l and is given in Section [XT] 

SectionjH applies the results of SectionjSjto a frequency-based probabilistic analysis of quan¬ 
tized linear systems, which is the main theme of the present report. A quantized linear system 
is defined in Section IATI as a dynamical system in Z” with the transition operator T ■= Ro L, 
where L G is a nonsingular matrix, and i? : —>■ Z” is a quantizer which commutes 

with the additive group of translations of Z” and such that is a Jordan measurable set. 

An example of a quantizer is the roundojf operator i?* with = [—1/2,1/2)"^. The quan¬ 

tized linear {R, L)-system is a model of the spatial discretization of a linear dynamical system 
in fixed-point arithmetic with no overflow. Section 14.2! defines an associated algebra Q of C- 
quasiperiodic subsets of Z"^, with the matrix C G formed from nonzero integer powers 

of L. By splitting C into two submatrices C~ and £+ formed from negative and positive integer 
powers of L, respectively, we define a backward algebra Q~ of £“-quasiperiodic subsets of Z” 
and a forward algebra Q+ of £+-quasiperiodic sets. It turns out that the associated algebra Q 
is invariant with respect to both the transition operator T and its set-valued inverse T~^, and 
that the backward and forward algebras Q~ and Q+ are invariant with respect to T and 
respectively. Section [431 shows that a mes„ 2 -almost every nonsingular matrix L G is iter¬ 
atively nonresonant in the sense that its positive and negative integer powers form a nonresonant 
matrix C := ^ Hence, for any such matrix L, the associated algebra Q (in¬ 

cluding its subalgebras Q~ and Q+) consists of frequency measurable subsets of the lattice IR. 
Moreover, it turns out that for an iteratively nonresonant matrix L, the corresponding transition 
operator T is not only Q+-measurable, but also preserves the frequency F on the forward alge¬ 
bra Q+, thus allowing the quadruple (Z”, Qf, F, T) to be regarded as a dynamical system with 
an invariant finitely additive probability measure F. Section l4Al is concerned with quantization 
errors which are defined as the compositions Ek := E o for k = 1,2, 3,..., where the 
first quantization error E maps a point a: G Z” to the vector E{x) := Lx — T{x) G i?“^(0). It 
is shown that, under the condition that the matrix L is iteratively nonresonant, the quantization 
errors Ei, E 2 , E^,... are mutually independent, uniformly distributed on the set and 
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are Q+-measurable in the sense that for any positive integer N and any Jordan measurable set 
G C the following set is an element of the forward algebra Q+ and its frequeney is 

mes NnG: 

{a: G Z” : {Ek{x))i^k^N ^ ^ ■ (1-1) 

These properties of the quantization errors, whieh extend V.V.Voevodin’s results of lIT^ . are 
elosely related to the property that the transition operator T is not only Q+-measurable and 
preserves the frequeney F on the forward algebra Q+, but is also mixing. Henee, the quadruple 
Q'*', F, T) ean be regarded as an ergodie dynamieal system. Using the property that the 
quantization errors are uniformly distributed over we define a supporting dynamieal 

system in with the affine transition operator %{x) := Lx — p, where p := 
is the mean vector of the uniform distribution. It is more convenient to study the deviation of 
trajectories of the quantized linear system from those of the supporting system whose phase 
portrait is qualitatively similar to that of the original linear system. Section 14.51 shows that 
the vector {T^{x) - T’^{x))i^k^N e R^", which describes the deviation of positive semi¬ 
trajectories of the quantized linear and the supporting systems during the first N steps of 
their evolution from a common starting point x G Z*^, is an affine function of the quantiza¬ 
tion errors Fi(x),..., E]\f{x), thus clarifying the role of the latter. Moreover, the sequence 
iPix) - driven by the quantization errors, is a homogeneous Markov chain on the 

probability space (Z"^, Q+, F). Therefore, in view of the results of Sections 1441 and |431 for any 
iteratively nonresonant matrix L, the events, which pertain to the deviation of positive semitra¬ 
jectories of the discretized system and the original linear system in a finite number of steps of 
their evolution, are representable as frequency measurable quasiperiodic subsets of the lattice 
in (11.11) which belong to the forward algebra Q+. 

In contrast to the original linear system, the nonlinear transition operator T is, in general, 
neither surjective nor injective. There are “holes” x G Z" in the lattice, which are not reach¬ 
able for T, that is, T~^{x) = 0, and there also exist points x G Z" for which the set T“^(x) 
consists of more than two elements. Therefore, the negative semitrajectory (T“*’(x))fcj.i of 
the quantized linear system is a set-valued sequence with values from the class of finite (pos¬ 
sibly, empty) jubsets of the lattice. Section 14.61 defines aj:ompensating system as the quan¬ 
tized linear (i?, L“^)-system with the transition operator T := R o where the quantizer 
R is given by R{u) := R{u + (J„ -f L~^)p). This allows the preimage T~^{x) to be repre¬ 
sented as the Minkowski sum T^{x) + 'L]y{x) of the singleton T^(x) and a finite (possibly, 
empty) set SAr(x) C Z"^. The set SAr(x) is completely determined by the quantization errors 
Ei{x),..., E]\f{x) for the compensating system which are defined similarly to the quantiza¬ 
tion errors of the quantized linear {R, L)-system. Under the assumption that the matrix L is 
iteratively nonresonant, the quantization errors of the compensating system are mutually in¬ 
dependent, uniformly distributed on the set i?“^(0) and Q“-measurable in the sense that, for 
any positive integer N and any Jordan measurable set H C (i?“^(0))^, the following set is an 
element of the backward algebra Q~ and its frequency is mesNnH: 

{x G Z" : {Ek{x))i^k^N G if } G Q~. (1.2) 

Due to these properties, the set-valued sequence {T~ '^{x)-T ^(x))fcs.i consists of the Minkow¬ 
ski sums of the corresponding elements of the sequence (T“^(x) — T^(x))fcj.i and the set-valued 
sequence {—T^k{x))k^i, with the last two sequences being homogeneous Markov chains on the 
probability space (Z", Q“,F). As a corollary, a martingale property [fT4l is established for 
the sequence (| det L\^jj^T~^{x))k^i which describes the cardinality structure of the set-valued 
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negative semitrajectory of the quantized linear system. In view of the results of Section 14.61 
under the assumption that the matrix L is iteratively nonresonant, any event, pertaining to the 
deviation of negative semitrajectories of the discretized and linear systems in a finite number of 
steps of their evolution, is representable as a frequency measurable quasiperiodic subset of the 
lattice in (11.21) which belongs to the backward algebra Q~. 

Therefore, Section |4] can be summarized as follows. The phase portrait of the original linear 
system with a nonsingular matrix L G is distorted when the system is replaced with the 
quantized linear {R, L)-system. This distortion can be described in terms of the deviations of 
positive and negative semitrajectories of the quantized linear system from those of the original 
linear system. For an iteratively nonresonant matrix L (such matrices are typical), these devia¬ 
tions can be studied in the framework of the probability space (Z”, Q, F), since all the events, 
related to the deviations in a finite number of steps of system evolution, are representable as 
frequency measurable quasiperiodic subsets of the lattice belonging to the associated algebra 

Q. 

Section [5] considers a particular class of neutral quantized linear (i?, L)-systems on the lat¬ 
tice Z" of even dimension n = 2r, with a matrix L being similar to an orthogonal matrix with a 
nondegenerate spectrum. This class of systems is specified in Section [5 t 1 The special structure 
of the matrix L is used in Section 15.21 in order to establish a functional central limit theorem 
for the deviation of positive semitrajectories of the quantized linear system and the supporting 
system. Section [531 applies the above results to quantitative analysis of the phase portrait of a 
two-dimensional neutral quantized linear system with which the problem on rounded-off pla¬ 
nar rotations is concerned IH [HI . We compute the frequencies of some events pertaining 
to the phase portrait of the system and compare these theoretical predictions with a numerical 
experiment on a moderately large fragment of the lattice. 

The results of this report (some of them were partially announced in [0) can be used for 
rigorous justification of ad hoc models [|5l for spatial discretizations of dynamical systems. The 
apparatus of frequency measurable quasiperiodic sets, proposed in the present study, is also 
applicable to quantitative analysis of aliasing structures [jTllHl. 


2 A probability structure on the integer lattice 


2.1 Averageable functions 

In what follows, IR denotes the integer lattice in the real Euclidean space of n-dimensional 
column-vectors with the standard inner product x^y and Euclidean norm |x| := \/x^x, where 
(•)^ is the transpose. Eet M denote the linear space of bounded functions / : Z”^ —)■ M with the 
uniform norm 

ll/ll := sup |/(x)|. 

Denoting the class of nonempty finite subsets of the lattice 7R by 


n := {P C Z” : 0 < #P < +CX)} , 


( 2 . 1 ) 


with the number of elements in a set, we define a functional W : M x IT —M which maps 
a function / G M to its average value over a set P G 11: 


wu.P) 


1 

w 


xeP 


( 2 . 2 ) 


6 











This quantity is also representable as the expeetation E/(,^), where is a random veetor with 
uniform probability distribution over the set P. For any vectors a := {ak)i!^k!in £ and 
N", with N the set of positive integers, let 

Pa,i := {a: := [xk)\<;ji^n G Z” : ^ a:fc < Ofc + 4 for all 1 ^ ^ n} (2.3) 

denote a discrete parallelepiped which consists of 11^=1 points of the lattice. For any given 
G N, we also define a class 

Vn := {Pa,i : a, i e IP, min 4 ^ ivl (2.4) 

( l^fcsgn J 

of sufficiently large parallelepipeds whose edge lengths are bounded below by iV. These sets 
form a decreasing sequence: Vi ^ V 2 ^ with the class 

p ;= ; N e N} (2.5) 

specifying a topological filter base on the set If in (12.11) . For any function / G M, consider the 
lower and upper limits of the average W{f,P) in (12.21) as the set P tends to infinity in the sense 
of (1231): 


A4/) := liminflP(/,P) := sup inf W{f,P), (2.6) 

P/'oo P&Vn 

A*{f) := limsuplfo(/,P) := inf sup W{f,P). (2.7) 

p/'oo P£Vn 

Definition 1 For a given function / : —)■ R, the limits (12.(51) and (12.71) are called the lower 
and upper average values of the function, respectively. If these limits coincide, their common 
value is written as A(/) and is called the average value of f, and the function is called aver- 
ageable. 

Therefore, the average value of an averageable function / is the limit A(/) := limp W (/, P). 
Constants are trivial examples of averageable functions. In particular, the identically unit func¬ 
tion 1 : Z" —)■ {1} is averageable, and its average value is A(l) = 1. Less trivial averageable 
functions will be studied in Section [3l 

Lemma 1 

(a) The functionals A*, A* : M —)■ R are concave and convex on the space M, are mono¬ 
tone with respect to the set M+ of nonnegative-valued bounded functions, are positively 
homogeneous and satisfy a Lipschitz condition, 

|A4/) - A4(7)|, |A4/) - A\g)\ ^ A4|/ - g\) ^ ||/ - gl (2.8) 

(b) the class of averageable functions M is a linear subspace of M; 

(c) the functional A ; Af ^ R (5 linear, positive with respect to the set Ai^ of 
nonnegative-valued averageable functions and satisfies |A(/)| ^ A*(|/|) ^ ||/||. 
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Proof. From (12.21) . it follows that for any given nonempty finite set P C Z", the funetional 
VF(-, P) : M —)■ M is linear, positive with respeet to the eone M+ and satisfies \W{f, P)\ ^ 
W{\f\,P) ^ ||/||. These properties imply the assertion of the lemma. ■ 

The average value funetional A is extended in a standard fashion to the eomplex linear spaee 
A4 + iA4 as A(/) := A(Re/) + iA(lm/). Furthermore, if / : Z” —)■ A is a map with values 
in an r-dimensional real linear spaee X, given by / := Ylk=i fA, where /^ : Z” — )■ R are 
averageable funetions and bi, ... ,br is a basis in X, then the eorresponding average value is 
defined to be A(/) ;= Ylk=i ^{fk)bk- This definition does not depend on a partieular ehoiee 
of the basis in X. 


2.2 Frequency measurable sets 

For any set A G let 


F*(A) A,(Xa) - hminf ^ , 

p/'oo jbP 

(2.9) 

F-(yl) := A*(Ix) = limsup * 

P/'oo 

(2.10) 

denote the lower and upper average values of the indicator function Xa{x) : = 

of the set. In (12.91) and (12.101). use is made of the equalitv 

J1 for X E A 
|o for X ^ A 


(2.11) 


for the relative fraetion of points of the set A in (a diserete parallelepiped) P, which follows 
from (12.21) . Note that 


0 ^ F*(A) = 1 - F*(Z” \ A) ^ F*(A) ^ 1. 

Definition 2 For a given set A C Z”, the quantities F* (A) and F* (A) in ( 12.91) and ( 12.701) are 
called the lower and upper frequencies of the set, respectively. If they coincide, then their com¬ 
mon value is denoted by F(A) and is called the frequency of A, and the set is called frequency 
measurable. 

Therefore, the frequency measurability of a set A is equivalent to the averagability of its indi¬ 
cator function X^. The frequency F(A) of a frequency measurable set A is the limit of the ratio 
(12.111) in the sense of (I2.5I) - (I2.7I) : 


F(A) : = 


lim 

P/'oo 


#(An^) 


That is, F(A) is the limit fraction of points of a given set A in unboundedly increasing rectan¬ 
gular fragments of the lattice. A straightforward example of a frequency measurable set is Z", 
with F(Z"^) = 1. Also note that any finite subset of Z" has zero frequency. A class of non¬ 
trivial frequency measurable sets, relevant to the context of spatial discretizations of dynamical 
systems, will be considered in Section [3l 
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Lemma 2 


(a) For any disjoint sets A, B G IF, their lower and upper frequencies satisfy the inequalities 

F,(A1Js) ^F,{A) + F,{B), F*(/Iljs) ^F*{A) + F*{B). 

Furthermore, F*(y4) ^ F*(i?) andF*{A) ^ F*{B) hold for any sets A G B G HA; 

(b) the class S of frequency measurable subsets oflA contains the lattice and is closed under 
the union of disjoint sets and complement of a set to the lattice; 

(c) the frequency functional F ; S —)■ [0,1] is a finitely additive probability measure, which 
satisfies F(Z”) = 1 and F(y4 [J i?) = F(74) + F(i?) for any disjoint sets A,Be S; 

(d) if the symmetric difference AAB of sets A, B G HA is frequency measurable with zero 
frequency F(74Ai?) = 0, then A is frequency measurable if and only if so is B, with 
F(y4) = F(S) in the case of frequency measurability. 

Proof. The assertion (a) of the lemma follows from the concavity, convexity, monotonicity 
and positive homogeneity of the functionals A* and A* (see the assertion (a) of Lemma [T]) and 
from the equality 

Aa\jb=Aa + 'Ab ( 2 . 12 ) 

for any disjoint sets A, B G HA. The assertions (b) and (c) of the lemma follow from the 
linearity of the space A1 and of the functional A (see the assertions (b) and (c) of Lemma [B, 
and from the equality (12.121) for disjoint sets. The assertion (d) of the lemma is established by 
using the inequalities 


|F,(A) - FfiB)\, |F*(A) - F*{B)\ ^ F*{AAB) 

which follow from (12.81) and from the identity Faab = \Aa —Ab\, thus completing the proof of 
the lemma. ■ 

The assertions (b) and (c) of Lemma [2] show that the triple (Z”, S, F) can be regarded as a 
probability space, with the average value functional A : A1 —)■ M playing the role of expec¬ 
tation. In view of Lemma [2td), this probability space is complete. However, this space does 
not satisfy the axiomatics of A.N.Kolmogorov [|T4ll since S is not a a-algebra and F is not 
a countably additive measure. Despite this pathology, frequency measurable sets and the fre¬ 
quency functional on them possess geometric properties which correspond to those of Lebesgue 
measurable sets and the Lebesgue measure in M". 

Theorem 1 For any set A G iS, any nonsingular matrix F G and any vector z G H"^, the 
set FA + z is also frequency measurable, with 
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Proof. We will first verify the equality (12.131) for the identity matrix F = In of order n. This 
is equivalent to that the translated set A + z inherits frequeney measurability from ^4 G 5, with 

F( 2 l + ^) = F( 2 l) (2.14) 


for any z ^ IF. Note that Pa/ + z = Pa+z,e for any a e and i G N”, where Pa/ is the 
diserete parallelepiped defined by (12.31) . Henee, eaeh of the elasses Vn in (12.41) is invariant 
under translations of the eonstituent parallelepipeds: 


Fn -f z {P z : P E. Pn} — Pn- 


In combination with the identities #((^4 + z) f) P) = #(^4 f)(P — z)) and — z) = #P, 
the translation invariance implies that 


inf 

P&Vn 


#((A+^)np) 

*p 


sup 

P&Vn 


#{{A + z)[]P) 


inf 

P&Pn 


*P 


sup 

P&Vn 


*p 


From the arbitrariness of G N in the last two equalities and from the assumption that y4 G 5, 
it follows that F *(74 + z) = F *(74 + 2 ;) = F( 74 ), which implies (12.141) . Therefore, it now 
remains to prove the property (12.131) for 2 : = 0. More precisely, we will show that the set PA is 
frequency measurable and 

S 

for any A G 5 and nonsingular P G To this end, for every u E N and a := {oik)i^k^n ^ 

IF, consider a discrete cube 


P^'^ := {x := {xk)i^k^n e : vau ^ Xk < i^iak + 1 ) for all 1 < A; ^ n} (2.16) 

which consists of i/” points. For any given u, the cubes Pa'* belong to the class P^ in (12.41) and 
form a partition of Z". Therefore, the sets A^, „ = A f) Pa\ a E I"', form a partition of the set 
A C Z"^, and for any nonempty finite set Pel"", 


F U A..C(FA)f]PcF U A.,a, (2.17) 


where 

$^(P) ;= {a G Z" : PP^^) C P} C ^^(P) := |a G Z" : (PP^^^^) p|P ^ 

In view of the assumption that det P 7 ^ 0, it follows from (12.171) that 

|/”#4„(P) inf AAM ^ #((F.4)f|P) < n“#>I>„(F) sup AAM. 

QeP^ HQ Q^-p^ HQ 


(2.18) 


We will now prove that the limits 


lim 

Ppoo 




lim 

pAoo 


#P 


1 

detP 


(2.19) 
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in the sense of (I2.5I) - (I2.7I) hold for any i/ G N. To this end, eonsider the Minkowski sums 
Pa,e '■= Pa/ + V and := Pi^^ + ^ of the diserete parallelepipeds Pa/, Pt'^ C Z” with 
the eube V := [—1/2,1/2)^^ C M"'. By using the matrix norm |||P||| := rnaxis;j,g„ \ fjk\, 

whieh is indueed by the Chebyshev norm in R”, with fjk denoting the entries of the matrix F, 
it follows that 


and 


1)11X1 +1)12) C y (FPM) 

(2.20) 

a&^u{Pa,i) 


U (FPM)CX,, + ((2P-1)11X1-1)1/ 

(2.21) 




where G-P = {mgR": u + H g G} denotes the Minkowski subtraction of sets G,H G R"^. 
The inclusions (12.201) and (12.211) imply that 


. i>''#'S>/p) ^ 1 ( _ (2-^-i)iiifiii + n 

PSV„ #P ^ |detf| V N ) 

pj„ #F ^IdetFlV N ) 


These inequalities lead to the limits in (12.191) . Now, from (12.181) and (12.191) . it follows that for 
any z/ G N, 


——— inf 

detP Q&v, 


#Q 


^ F*(P2l) ^ F*(P2l) ^ 


det PI 


sup 


moQ) 


and hence. 


F,(21) 
I detP 


^ F*(P2l) ^ F*(P2l) ^ 


F*(2l) 

|detP|‘ 


The last inequalities and the assumption A G 5 imply (12.151) . thereby completing the proof of 
the theorem. ■ 


A geometric interpretation of (12.131) is that the frequency functional F is translation invari¬ 
ant, and a linear transformation of a frequency measurable set A with an integer matrix P, 
satisfying | det P| > 1, leads to a “sparser” set FA. 


2.3 Distributed maps 

Suppose X is a metric space equipped with the Borel cr-algebra jB, and let P be a countably 
additive probability measure on the measurable space (X, B). Recall that a map g : X ^ Y 
with values in another metric space Y is said to be P-continuous [[U if the discontinuity set 
of g has zero P-measure. Accordingly, a set P G P is called P-continuous if its boundary 
dB (which is the set of discontinuity points of the indicator function Ib of the set B) satisfies 
D{dB) = 0. 

Definition 3 Suppose g ■. IF ^ X is a map with values in a metric space X, and D is a 
countably additive probability measure on (X, B). The map g is said to be D-distributed if the 
preimage 

g-\B) := {x G Z- : g{x) G B} 
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of any D-continuous set B ^ B under the map g is frequency measurable and its frequency 
satisfies 

F{g-\B)) = D{B). (2.22) 


In this case, the algebra 


Sg := {g-\B) : Be B, D{dB) = O} , 

which consists of frequency measurable sets, is referred to as the algebra generated by the map 

9- 

Lemma 3 Suppose g ■. IP ^ X is a map with values in a metric space X, and let D be a 
countably additive probability measure on (X, B). Then the following properties are equiva¬ 
lent: 

(a) for any bounded continuous function f : X ^ M., the composition / o ^ —)■ R 

averageable, with 

Mf°9)= [ fix)D{dxy, (2.23) 

Jx 

(b) for any bounded D-continuous function / : X —)■ R, the composition f og is averageable 
and ( 12.231) holds; 

(c) the map g is D-distributed. 

Proof. For any nonempty finite set P C Z**, eonsider a eountably additive probability measure 
Up on (Z”, 2^") defined by 


Up{A)-.= W{Xa,P) 


#P 


where use is made of (12.21) . Any such set P and any map g -.IP ^ X generate a countably ad¬ 
ditive probability measure Upo g~'^ on (X, B). The expectation of a Borel measurable function 
/ : X —>■ R, interpreted as a random variable on the probability space (X, B, Up o g~^), takes 
the form 


f{x)(Upog-fi[dx) = WUog,P). 


Jx 

In particular, application of the equality (12.241) to the indicator function / 
yields 


(PpO^-l)(P) 


#{g-\B)C]P) 

#P 


(2.24) 

Xp of a set P c X 

(2.25) 


Now, we assume the map g to be fixed but otherwise arbitrary, and, similarly to (12.41) and (12.51) . 
define a topological filter base 


D := {Pjv : N e N}, Vm := {Up o g-^ , Be Vn}, (2.26) 


on the class of countably additive probability measures on (X, B). The existence of a limit D in 
the sense of (12.261) in the topology of weak convergence |[I1 of probability measures on (X, B) 
means that 
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(a’) for any bounded continuous function / : X —)■ M, 


lim [ f{x){Upog ^)(da:) = [ f{x)D{dx). (2.27) 

Jx Jx 

By the well-known criteria for the weak convergence, the last property is equivalent to each of 
the following ones: 

(b’) the convergence (12.271) holds for any bounded D-continuous function / : X —)■ M; 

(c’) for any D-continuous set B E B, 


hm {Upog-^){B) = D{B). 

P/'oo 


(2.28) 


The properties (a’)-(c’) are equivalent to the corresponding properties (a)-(c) stated in the 
lemma. Indeed, in view of (12.241) . the left-hand side of (12.271) is A(/ o g), and by (12.251) . 
the left-hand side of (12.281) is F{g~^{B)). Hence, the equivalence of the properties (a’)-(c’) im¬ 
plies that (a)-(c) are also equivalent to each other, thereby completing the proof of the lemma. 

■ 

From the proof of Lemma[3l it follows that for any map g -.1/^ ^ X with values in a metric 
space X, there exists at most one countably additive probability measure D on (X, B) such that 
g is D-distributed. 

Lemma 4 Suppose Xi and X 2 are two metric spaces with Borel a-algebras Bi and B 2 , respec¬ 
tively. Also, let gi ■. IP ^ Xi and g 2 '■ Xi ^ X 2 be two maps which satisfy the following 
conditions: 

(a) gi is Di-distributed, with Di a countably additive probability measure on (Xi, Bi); 

(b) g 2 is Dpcontinuous. 

Then the composition of the maps g 2 o gi '■ IT ^ X 2 is D 2 -distributed, with the probability 
measure D 2 on (X 2 , B 2 ) given by 

D 2 {B) = D,{gf\B)), B e B 2 , (2.29) 

and the generated algebras satisfy the inclusion 

<Sg20gi C Sgj^. (2.30) 

Proof. Let 99 : X 2 — )■ M be a bounded continuous function. Then, by the assumption (b) of 
the lemma, the function 

/:= 9905(2 : XiM ( 2 . 31 ) 

is bounded and Di-continuous. Hence, by the assumption (a) of the lemma and by the criterion 
(b) of Lemma [3l the function / o 5 ]^ : Z" —>■ M is averageable and 

Mf°9i)= [ f{x)Di{dx)= [ (p{y)D 2 {dy), (2.32) 

JXi JX2 
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where D 2 is the probability measure given by (12.291) . Since the function cp is otherwise arbitrary, 
then, in view of the criterion (a) of Lemma [H it follows from (12.311) and (12.321) that the map 
92 o 91 is i 92 -distributed. It now remains to prove the inclusion (12.301) . From the assumption 
(b) of the lemma and from (12.291) . it follows that D 2 {dB) = Di{d 92 ^{B)) for any B G B 2 , 
and hence, 92 ^{B) is a i9i-continuous Borel subset of Xi for any i 92 -continuous set B G B 2 . 
Therefore, 

{9^\92\B)) : Be B2, D2{dB) = 0} C {9^\B) : B e B,, D,{dB) = 0}. 

Since the left-hand side of this inclusion is the algebra Sg^ogi, whilst the right-hand side is iS^j, 
the inclusion (12.301) follows. ■ 


3 Quasiperiodic objects on the integer lattice 

3.1 Cells 

In what follows, we will use an extension of the notion of a space-filling polytope [|3]| given 
below. 

Definition 4 A set V C is called a cell if its translations V + z, considered for all z G Z'", 
form a partition of MB'. 

An example of a cell in R'" is provided by the cube [0,1)'". Recall that any unimodular matrix 
(that is, a square integer matrix with determinant ± 1 ) describes a linear bijection of the integer 
lattice. 


Lemma 5 

(a) For any unimodular matrix F G any u G R"^ and any cell V C R"^, the set 

FV + uis a cell in R*”; 

(b) a set V C R'" is a cell if and only if there exists a partition 

{fl, : z G IB} (3.1) 

0 / [0,1)"^ satisfying 

y=lj(n, + ;^). (3.2) 

Moreover, such a partition is unique; 

(c) any Lebesgue measurable cell V C R"^ has a unit measure, mesml^ = 1; 

(d) for any cells If C R”^^ and If C R"^^, the set If x If is a cell in R'^i+^^a^ 
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Proof. To establish the assertion (a), note that any unimodular matrix F G determines 

a linear bijection of Z™. Henee, for any eell V C any u E M'” and any z G IF \ {0}, the 
set G := FV + u satisfies 


G + Z^ = u + F{V + F-^Z^) = R'", 

Gf]iG + z)=u + F(vf]iV + F-^z)'^ =0 

which means that G is also a cell. In order to prove the assertion (b), we associate with a given 
but otherwise arbitrary set V C R'" the sets 

ZGZ™, (3.3) 

in terms of which the set V is representable by (13.21) . that is, V = + ^)- Indeed, for 

any 2 ; G Z”", 

fl° + ;2=([0,l)”^ + z)flf", (3.4) 

and hence, 

zeZ"* VzeZ" / 

On the other hand, the sets 0° in (13.31) partition the cube [0,1)'" if and only if 1/ is a cell in 
R'". This proves the first part of the assertion (b). We will now prove that the sets 0° in (13.31) 
provide a unique partition of [0,1)”^ such that a given cell V C R™ is represented by (13.21) . To 
this end, let (13.11) describe an arbitrary partition of the cube [0,1)™ satisfying (13.21) . Then for 

any 2 ; G Z™, 

n, + z = {[o,ir + z)f]v, 

which, in view of (13.31) . implies that = 0°, thus proving the uniqueness of the partition. The 
assertion (c) follows from the relations 

meSrnV^ = ^ mesm(0° + z) = ^ meSmO° = mesm[0,1)™ = 1 

2eZ"* 

which hold for any Lebesgue measurable cell V C R"* and are based on (13.31) . (13.41) and the 
translation invariance of the Lebesgue measure. The assertion (d) of the lemma follows from 
the identities ^ •^m 2 ]gmi+m 2 _ ^ ]gm 2 _ g 

Lemma [5] shows that there exist more complicated cells in R'” than the cube [0,1)”^. This is 
illustrated by Fig. [T] which provides an example of such a cell in R^. 

Lemma 6 

(a) Let f : R™ -E- M be a locally integrable function, unit periodic with respect to its m 
variables: 

f{u + z) = f{u) forallM G R”^, ^ G Z™. (3.5) 

Then fy f(u)du = f(u)dufor any Lebesgue measurable cell V C R'"; 

(b) any translation invariant set G C R"^ (under the group of translations ofZF in the sense 
that G + Z'" = G) is representable as G = (V^ f) G) + Z'" where V is an arbitrary cell 
in R™; 
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Figure 1: An example of a eell in whieh is obtained by eutting the square [0,1)^ into pieees 
and translating them by two-dimensional integer veetors. 

(c) Let G C R™ be a translation invariant and Lebesgue locally measurable set. Then, for 
any Lebesgue measurable cell V C R'", 

mes^ n ^ n • 

Proof. The assertion (a) of the lemma is proved by the following equalities for any unit peri- 
odie and loeally integrable funetion / : R™ —)■ R: 

[ f{u)du = ^ [ f{u)du = ^ [ f{u)du = [ f{u)du, 

whieh employ a partition (13.11) from Lemma [2b) for a Lebesgue measurable eell V C R™ 
and the translation invarianee of the Lebesgue measure. The assertion (b) of the lemma is 
established by the relations 

(l/f|G')+Z™= IJ (^{V + z)f]{G + z)^ ={V + Z^)f]G = G 

whieh hold for any eell V C R'" and any translation invariant set G C R”^. The assertion (e) of 
the lemma follows from the assertion (a) and from the property that the indieator funetion Iq 
of any translation invariant set G C R™ is unit periodie with respect to its m variables. ■ 

3.2 Quasiperiodic sets 

With any positive integer m E N, any set G C R”^ and any matrix A E we associate a 

set 

Qm{G, A) := {x e Z" : Ax E G + Z™} (3.6) 

The following lemma provides a useful formalism for set theoretic operations with such subsets 
of the lattice Z"^ which will play an important role in the subsequent sections. 

Lemma 7 The subsets o/Z", defined by ( I3.<5D . possess the following properties: 
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(a) for any m eN, any A G any G C M™' and any cell V C 


QUG, A) = (V + Z™), a) ; 


(b) for any m eN, any A E any cell V C M'" and any G, H C V, 

Qm{V,A)=Z^, 

Z^\Q^{G,A) = Qm{V\G,A), 
QUO, A) f| QUH, A) = Qm{Gf] H, a) , 
QUO, A) IJ Qm{H, A) = Qm{G[j H, a) ; 

(c) for any m G N, G C and A E 


Qn+m ( G, 


= Q^{H,A), 


where 


H ■= {vE 


there exists u E'L'^ such that 


eG ; 


(d) for any m eN, any G C R™, any A E and any unimodular F E 

Q™(FG,FA) = g^(G,A); 

(e) for any mi, m2 G N, Gi C R^^h Ga C R™h Ai G R”"!''", Aa G 


Qm\-\-m.2 ( ^ Ga, 


Ai 

Aa 


Qmi (Gi, Ai) Qm2{G2-i Aa); 


(3.7) 


(3.8) 


(f) for any m G M, A G G C R”^ and z G Z”', 

Qm,{G, A) + ^ = Qm{G + Az, A). 


Proof. Since the set G + Z'" is translation invariant, then, in view of Lemma [6l)b), it ean be 
represented as 

G + Z^ = (v p|(G + Z™)) + Z”^ 

for any eell V C R™. In view of (13.61) . this implies the assertion (a) of Lemma|7l The assertion 
(b) follows from the property immediately below. For any given eell V C R™, the map K : 
2 V 2 *'", defined by K{G) = G + Z'", satisfies K(V) = R™ and preserves the set theoretie 
operations in the sense that 


R”^\iF(G) = K{V\G), 

K{G)f]K{H) = K{Gf]H), 

k(g)\Jk{h) = k(g\Jh) 


17 










belongs 

y e 

G G. The last three relations hold if and 


for any G,H C V. In order to prove the assertion (c), we note that a point x G 
to the set on the left-hand side of (13.71) if and only if there exist m G M”, n G S 

11 

z G Z”* satisfying x = u + y,Kx = v-\-z and ^ 

only if X belongs to the set on the right-hand side of (13.71) . with the set H given by (13.81) . The 
assertion (d) follows from the property that {FG) + TF = F{G ^ Z”*) for any unimodular 
matrix F G and any set G C M™. The assertion (e) of the lemma follows from the 

equality ((S*! x G 2 ) + + Z™i) x {G 2 + for any sets Gi C and 

G 2 C Finally, the assertion (f) is proved by noting that for any z G a point x E ZP 
satisfies Kx E G + ZT^ if and only if A(a: P z) E G + Kz + Z"*. ■ 


By Lemma |7] (a), we can restrict ourselves, without loss of generality, to considering the 
sets Qm{G^ A) in (13.61) only for G C V, where L is a fixed but otherwise arbitrary cell in MF. 
For any matrix A G we define the class of sets 


Qm{^) '■= {Qm{G^ A) : (7 C [0,1)”^ is Jordan measurable} . (3.9) 

From Lemma |7}b) and the property that Jordan measurable subsets of [0,1)”* form an algebra, 
it follows that Qm(A) in (13.91) is an algebra of subsets of Z"^. 


Definition 5 The algebra Qm{^), defined by ( 13.91) for given m E N and A G is called 

the algebra of A-quasiperiodic subsets o/Z". 


Lemma 8 The algebra (13.91) possesses the following properties: 


(a) for any unimodular matrix F E Qm{FA) = Qm{A); 

(b) if Ai E is a submatrix of A 2 E then Qmi(Ai) C Qm 2 {^ 2 )i 

(c) for any m G N and any A E R™'^"', 


Q 


n+m 



Qm{A)- 


(3.10) 


(d) the algebra Qm{A) is translation invariant in the sense that A E Qm(A) implies A + z E 
Qm{A) for any z G Z”^. 


Proof. In order to prove the assertion (a), we note that, in view of the assertions (a) and (d) of 
Lemma |71 

Q^{G,FA) = Q^{K{G),A) 

holds for any unimodular matrix F G and any set G C [0,1)'". Here, the set K{G) C 

[ 0 , 1 )*" is given by 

K{G) = [0,1)”^ + Z^). 

The map Tf is a bijection on the algebra of Jordan measurable subsets of [0,1)'", whence the 
assertion (a) of the lemma follows. In order to establish the assertion (b), it suffices to consider 
the case where 
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with A 3 G ]R(”^ 2 -mi)xn Indeed, this structure can always be achieved by permuting the rows 
of the matrices Ai and A 2 , which, in view of the assertion (a) of the lemma, does not affect the 
corresponding algebras of quasiperiodic sets Qmi(Ai) and Qm 2 {^ 2 )- By using the assertions 
(b) and (e) of Lemma |71 it follows that any A 1 -quasiperiodic set Qmi{G,ki) (with a Jordan 
measurable set G C [0,is representable as 

Ai) = X [0, ir™, A 2 ). 


The right-hand side of this equality is a A 2 -quasiperiodic subset of thus establishing the 
assertion (b) of the lemma. In order to prove the assertion (c), we note that, since A G is 


a submatrix of 


, then, in view of the above assertion (b), the algebra Qm(A) is a subalgebra 
of the algebra on the left-hand side of (13.101) . On the other hand, by applying Lemma |7tc), it 
-quasiperiodic set is representable as a A-quasiperiodic set. Hence, these 


follows that any 


two algebras also satisfy the opposite inclusion and the equality (13.101) . The assertion (d) of the 
lemma follows from Lemma IT^f)- • 

For what follows, we need an extension of the algebra of quasiperiodic sets (13.91) to infinite¬ 
dimensional matrices A. More precisely, for a given matrix 


C := 


(A, 


‘jkJjGX, l^k^n 


we define 


Qoo(^) ■= {Qm(G, A) : G E [0,1)™' is Jordan measurable, 
A G is a submatrix of £, m G N} 

= U U S„(A) 

m^l A is an (mxn)—submatrix of £ 

which is also an algebra of subsets of the lattice Z"'. 


(3.11) 


(3.12) 


Definition 6 The algebra Qoq{G), defined by (13.721) for a given matrix C G is called an 

algebra of C-quasiperiodic sets. 


3.3 Frequency measurability of quasiperiodic sets 

With any a; G we associate an elementary trigonometric polynomial —)■ C defined 

by 

T^(a;) := e^™ ", (3.13) 

where i := is the imaginary unit. 

Lemma 9 For any u G R", the function in (13.731) is averageable and its average value is 

= Xgn (cj). 
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Proof. If a; G then T^(a;) = 1 for all x G Z", and hence, A{T^) = 1 = Xzn(a;) holds for 
such vectors u. Now, suppose u := {ujk)i^kiin ^ X", that is, the indices of noninteger entries 
of u form a nonempty set 

/C := {1 ^ ^ n : ccfc ^ Z} . (3.14) 

The function in (13.131) has the following average value (12.21) over the discrete parallelepiped 
Pa,i in (1231): 


W(T^, Pa,i) = exp 2m ^ Wfc I Ufc + 


(3.15) 


k=l ^ ' J k=l 

where the function -0 : N x M —)■ [—1,1] is defined by 

((—for u G Z 

^{u,v) = I 

for u^Z 

V usin(7rf;) ^ 

From (13.151) and (13.161) . it follows that for any iV G M, 

sup \wmp )\«^ n I. / M . 

p&Vm I sm(7ra;fc) | 


(3.16) 


(3.17) 


where is the class of sufficiently large parallelepipeds given by (12.41) . Since the set /C in 
(13.141) is not empty (and hence, #/C > 0), then the right-hand side of (13.171) converges to zero 
as N + 00 . This convergence implies that, for any vector u G M” \ Z”^, the function in 
(13.131) is averageable, with A{TiS) = 0 = Izn{uj), thus completing the proof of the lemma. ■ 


Definition 7 A matrix A G is said to be nonresonant if the rows of the matrix 

]^(n+m)xn linearly independent over the field of rationals. 


G 


Note that nonresonant matrices do exist. Moreover, resonant matrices A G R™^” (which are 
not nonresonant) form a set of mn-dimensional Lebesgue measure zero. Indeed, since the 
nonresonance property is equivalent to 

A^m^Z" forallMGZ™\{0}, (3.18) 

then the set of all resonant matrices A G R"^^" can be represented as a countable union 

I I r 

IJ u,v 

of (m — l)n-dimensional affine subspaces := {A G R'"^" ; A^u = v}, each of which has 
zero m?7,-dimensional Lebesgue measure. 


Theorem 2 Suppose A G R'"^” is a nonresonant matrix, and a function f : R"^ —)■ R satisfies 
the conditions 


(a) / is unit periodic with respect to its m variables; 


(b) / is mesm-continuous; 


(c) / is bounded. 

Then the function / o A : Z" —>■ R w averageable, and its average value is computed as 

A{foA)= [ f{u)du. (3.19) 
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Proof. We will follow a standard scheme which is known as the method of trigonometric sums 
in number theory [|T^ and as Weyl’s equidistribution criterion in the theory of weak convergence 
of measures [[I|. More precisely, the proof will be carried out in three steps: we will establish 
(13.191) for trigonometric polynomials /, then for continuous functions / and finally, for arbitrary 
functions / satisfying the assumptions (a)-(c) of the theorem. Throughout the proof, A G 
is a nonresonant matrix in the sense of (13.181) . 

Step 1. Suppose / : M"* —)■ C is a trigonometric polynomial, that is, a linear combination 
of functions from (13.131) : 

f-.= Y,Cu:T^- (3-20) 

Here, c^j are complex coefficients, and is a finite subset of 7/^, which, without loss of gener¬ 
ality, is assumed to contain the zero vector. Then 


Co = / fiu)<lu. (3.21) 

-'[ 0 , 1 ]’" 

The composition / o A of the function / with the linear map specified by the matrix A is also a 
trigonometric polynomial: 

/ O A = ^CojTat^. 

By using Lemmas [H and |9l it follows that the function / o A is averageable with the average 
value 

A(/oA) = 5 ^c^A(Tat,) = 5^c^Xz.(ATa;). (3.22) 

Under the nonresonance condition (13.181) for A, the relation (13.221) reduces to A(/ o A) = cq, 
which is equivalent to (13.191) in view of (13.211) . 

Step 2. Suppose / : —)■ R is a unit periodic continuous function. By the Weierstrass 

approximation theorem, for any £ > 0, there exists a trigonometric polynomial : R'" —)■ C, 
defined by (13.201) . such that 

max \f{u)- fe{u)\^e. 

Therefore, by using the Lipschitz continuity of the upper average value functional (see Lemma[TJa)) 
and Step 1, 


A*(/o A) ^ A(/£ o A) + e = [ fe{u)du + e^ f f{u)du + 2e. (3.23) 

-'[ 0 , 1 ]"’ -'[ 0 , 1 ]"’ 


In view of the arbitrariness of £ > 0. it follows from (13.231) that 


A*{foA)^ j 

f f(u)du. 

[O,!]"* 

(3.24) 

A similar reasoning leads to 



A*(/ ° A) ^ j 

1 f(u)du. 

[O,!]"* 

(3.25) 


The inequalities (13.241) and (13.251) imply the averagability of the function / o A, with the average 
value given by (13.191) . 


21 


























Step 3. Now, let / : R"* — )■ R be an arbitrary funetion satisfying the eonditions (a)-(e) of 
the theorem. Consider the functions : R"^ —>■ R defined by 

f~{u) = min ( /(m), liminf f{v)) , = max ( f{u), limsup f{v) ) . 

\ V / \ / 


Both functions / and /+ are unit periodic in their m variables, and are lower and upper semi- 
continuous, respectively. Also, they satisfy the inequalities 


/ {u) ^ f{u) ^ /+(m) 


which turn into equalities for mesm-almost all u G [0,1)'" in view of the mes^-continuity of 
the function /. Hence, there exists a decreasing sequence of unit periodic continuous functions 
: R"^ —)■ R which converge to the function point-wise in the cube [0,1)'” and hence, to 

the function / almost everywhere in [0,1)™ as A; —)■ -|-oo. Therefore, by using Step 2, it follows 

that 

A*(/oA)^ f fk{u)du. (3.26) 

Application of the Lebesgue dominated convergence theorem to the right-hand side of (13.261) 
leads to the upper bound (13.241) . The lower bound (13.251) is verified similarly. As before, (13.241) 
and (13.251) imply the avaragability of the function / o A together with (13.191) . thus completing 
the proof of the theorem. ■ 


Theorem 3 For any nonresonant matrix A G and any Jordan measurable set G C 

R"*, the A-quasiperiodic set Qm{G,A) in (13.(51) is frequency measurable, and its frequency is 
computed as 

F(Q^(G, A)) = mes^ {v f|(G + IT)) , (3.27) 

where V C R'" is an arbitrary Lebesgue measurable cell. 

Proof. The Jordan measurability of the set G C R™ ensures that the indicator function 
/ := Xg+i.^ satisfies the conditions (a)-(c) of Theorem Hence, if the matrix A G R'"^" 
is nonresonant, then Theorem [2] implies that the function 

/ o A = Xg+z™ o a (3.28) 


is averageable, and 

A(/ o A) = / f{u)du = mes™ ([0,1)™ 0(67 + Z^)) . (3.29) 

Therefore, since the right-hand side of (13.281) is the indicator function of the set QmiG, A) in 
(13.61) . then this set is indeed frequency measurable, and, in view of (13.291) . its frequency is given 
by (13.271) for any Lebesgue measurable cell V C R'", with the last property following from 
Lemma [^c). ■ 

According to Theorem[3l if A G R™^"^ is a nonresonant matrix, then all the A-quasiperiodic 
sets, which form the algebra Qm{A) in (13.91) . are frequency measurable subsets of the lattice 
This result can be generalised to infinite-dimensional matrices C in (13.111) by appropriately 
extending Definition |71 
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Definition 8 A matrix C G is said to be nonresonant if all its submatrices A G 

are nonresonant in the sense of Definition^for any m G N. 

The existence of nonresonant matrices C G is established similarly to that in the case of 

finite-dimensional matrices. 

Theorem 4 Suppose a matrix C G is nonresonant. Then the algebra Qqo{C,) of C- 

quasiperiodic sets in i lAi2D has the following properties 

(a) Qoo('C) consists of frequency measurable subsets of IT; 

(b) for any N G ^ and any collection of N pairwise nonoverlapping submatrices Ck G 

of the matrix C, the corresponding algebras Qoo{C,k) cire mutually independent in the 
sense that 



for any Ak G Qoo(A)- 

Proof. The assertion (a) follows from Theorem [3] and from the property that for any A G 
Qoo(-^)> there exists a finite-dimensional submatrix A G R™^” of the matrix £ such that A G 
Qm(A). In order to prove the assertion (b), we note that for any sets Ak G Qoo(£fc) described 
in the theorem, there exist finite-dimensional submatrices A^ G of the corresponding 

matrices Ck (and hence, Ai,..., Aat are also non-overlapping submatrices of the matrix £), 
and Jordan measurable sets Gk C [0, such that Ak = Qm^i^ki^k)- Therefore, repeated 
application of Lemma |3e) yields 


N 

( 


’Ai' 

\ 


Gi X . 

• X Gn, 



k=l 



An 

/ 


(3.30) 


Since the matrix C is nonresonant, then the right-hand side of (13.301) is a quasiperiodic set (asso¬ 
ciated with a nonresonant matrix) whose frequency measurability is guaranteed by Theorem [3l 
with 

( N \ N N 

\ = JJ meSm^,Gk = JJ F{Ak), 

k=l J k=l k=l 

thus completing the proof of Theorem |4l ■ 


3.4 Distributed quasiperiodic maps 

Theorem 5 Suppose a matrix A G R™^"^ is nonresonant in the sense of Definition^ Also, let 
g : R™ X be a map with values in a metric space X such that 

(a) g is unit periodic with respect to its m variables; 

(b) g is mesm-continuous. 
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Then the composition goK ■.7T' ^ X is D-distributed, with D a probability measure on {X, B) 
given by 

D{B) := mes^ (vf]g-\B)'^ , (3.31) 

where V C is an arbitrary Lebesgue measurable cell. Furthermore, the algebra SgoA, 
generated by the map g o A, consists of A-quasiperiodic sets: 

SgoA C Qm{A). (3.32) 

Proof. For what follows, let niodm : K"* [0,1)"* denote the map whieh sends a veetor 

u := (Mfc)Kfe^m to the veetor 


modm(M) := (§Mfc}}) 

where {{■ J denotes the fraetional part of a number. The map mod^ is unit periodie with respeet 
to its m variables and identieally maps the eube [0,1)”* onto itself. We will now prove that if 
A G is a nonresonant matrix, then the eomposition mod^ o A : Z'” —)■ [0,1)*" is mesm- 
distributed. To this end, note that, for any bounded eontinuous funetion p : [0,1)™ —)■ R, the 
funetion f := p o modm satisfies the eonditions (a)-(e) of Theorem [2] and henee, the funetion 
/ o A is averageable, with 

A{p o modrn o A) = / f{u)du = / p{u)du. (3.33) 

i[0,l]’" -'[0,1]™ 

Sinee the funetion p is otherwise arbitrary, then, in view of the eriterion (a) of Lemma [3l it 
follows from (13.331) that the map modm o A : Z” —)■ [0,1)*" is mesm-distributed. Furthermore, 
applieation of Definition [3] yields 


‘^modrnoA — Qm{,A). (3.34) 

Now, suppose a map g : R'" —)■ X with values in a metrie spaee X satisfies the assumptions (a) 
and (b) of the theorem. In partieular, the unit periodieity of g implies that 

g o A = g o modrn o A. (3.35) 

Sinee the map g is mesm-continuous, and mod^ o A has been proved above to be mesm- 
distributed, then, in view of LemmaS the map o A is ^-distributed, where D is a probability 
measure on {X, B) given by 


D{B) =mesm ([0,1)™ f| ^"'(5)) . (3.36) 

Due to the unit periodieity of g, the preimage g~^{B) of any set i? C X is a translation invariant 
subset of R'". Henee, in view of the assertion (e) of Lemmathe right-hand side of (13.361) eo- 
ineides with that of (13.311) for any Lebesgue measurable eell V C R'". Moreover, by Lemma S 
it follows from (13.351) that SgoA C iSinod,„oA- The latter inelusion, eombined with (13.341) . implies 
(13.321) . thus eompleting the proof of the theorem. ■ 
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4 Quantized linear systems: general case 

4.1 Definition of a quantized linear system 

Suppose i? : —>■ Z" is a map which commutes with the additive group of translations of the 

lattice Z"^: 

R{u + z) = R{u) + z for all u G M”, z G Z”. (4.1) 

Such a map R is completely specified by the set which is a cell in in the sense of 

Definition m Indeed, in view of (14.11) . the preimages 

R-^{z) := {m G : R{u) = z} = R-\Q) + z, 

considered for all z G Z"^, form a partition of Moreover, for a given set 1/ C M", there 
exists a unique map i? : M" —)■ Z"^, satisfying (14.11) with = V, if and only if 1/ is a cell 

in M". 

Definition 9 A map i? : M" — )■ Z” is called a quantizer if it satisfies ( I4.il) and the set i?“^(0) is 
Jordan measurable. 


There is a one-to-one correspondence between quantizers and Jordan measurable cells in 
An example of a quantizer is as follows. 


Definition 10 The quantizer i?* : —)■ IT, with ^(0) = [—1/2,1/2)” is called the round¬ 

off quantizer. 


The roundoff quantizer ii* maps a vector u := {uk)i^kiin 
lattice Z” given by 


Rfiu) 



G R” to a nearest node of the 


with [-J the floor function. This is an idealised model of discretization in fixed-point arithmetic 
(with no overflow taken into account). In the framework of this model, there are more compli¬ 
cated quantizers R whose cells i?“^(0) are different from the cube [—1/2,1/2)”, as discussed 
in Section [3Tj 


Definition 11 Suppose R : R” —)■ Z” is a quantizer and L G R”^” is a nonsingular matrix. 
The dynamical system in the state space Z”, with the transition operator T -.JT given by 

T:=RoL, (4.2) 

is called a quantized linear {R, L)-system. 

The map T, defined by (14.21) . provides a model for the spatially discretized dynamical system 
which arises in simulating a linear system in R” (specified by the matrix L) on a computer with 
fixed-point arithmetic. 
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4.2 Associated, backward and forward algebras 

The following lemma shows that the elass of quasiperiodie subsets of the lattiee, introdueed 
in Seetion 13.21 is elosed under the dynamies of the quantized linear system in the sense that 
quasiperiodie sets are transformed to quasiperiodie sets. 

Lemma 10 For any m eN, any set G C and any matrix A G the quasiperiodie set 

Qm{G, A) in t l3.(5D is transformed by the transition operator T in d4.2D and its set-valued inverse 
T~^ as 


T-\Qm{G,A))=Qn+n. 


T{Qm{G, K)) — Qm+n 


({0} X G + 

'In 

A 

R-\0), 

'In 

A 

(g X {0} - 

'A' 

L-^R-\0), 




(4.3) 

(4.4) 


Proof. The definition (13.61) implies that a point x E IF belongs to the set on the right-hand 
side of (14.31) if and only if there exist 

ueG, veR-\ 0), yEl^, zEl^ (4.5) 


satisfying 


Lx = 


+ 


V + 


In terms of the eorresponding subveetors, the last equality is equivalent to 


Lx = V + y, A{Lx — v) = Ay = u -\- z. (4.6) 

On the other hand, in view of (14.51) . the leftmost equality in (14.61) is equivalent to T(x) = y, 
whereas the rightmost equality in (14.61) is equivalent toy E Qm (G*, A). Therefore, the fulfillment 
of (14.51) and (14.61) is equivalent to T(x) E Qm{G^A), that is, x E T~^{Qm{G^A)). Sinee the 
point X in the above eonsiderations was arbitrary, the set on the right-hand side of (14.31) indeed 
eoineides with T~^{QmiG, A)). Although the proof of the representation (14.41) is similar, we 
will provide it for eompleteness of exposition. A point x E I"' belongs to the set on the right- 
hand side of (14.41) if and only if there exist 

ueG, V e R~^{0), yEl"-, z eT^ (4.7) 

satisfying 


'A 


u 


'A' 

T ^ 

z 


L ^x = 


— 


L + 


A 


0 


A. 


y. 


The last relation is equivalent to 

Ly = v-\-x, A{L~^xL~^v) = Ay = u z. (4.8) 

On the other hand, in view of (14.71) . the leftmost equality in (14.81) is equivalent to Tiy) = x, 
whilst the rightmost equality in (14.81) is equivalent to y E Qm{G, A). Therefore, the fulfillment 
of (14.71) (14.81) is equivalent to x G T{Qm{G, L)). Hence, by the arbitrariness of the point x, it 
follows that the set on the right-hand side of (14.41) is equal to T{Qm{G, A)), thus completing 
the proof of the lemma. ■ 
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Lemma [TOl shows that the action of the transition operator T of the quantized linear system 
or its inverse T~^ on a quasiperiodic set QmiG, A) modifies the matrix A, thus leading to a set 
with a qualitatively different quasiperiodicity pattern. However, an algebraic closedness can be 
achieved here by restricting A G to submatrices of an infinite matrix C G which 

satisfies the property 


A £ 


In 

A 


L C and 


A 

In 


L-^ C, 


(4.9) 


where A <gC B signifies “A is a submatrix of B” for matrices with the same number of columns 
n. In particular, (14.91) implies that such a matrix C must contain L and L~^ as submatrices. 
Hence, by induction, the minimal matrix £, which satisfies (14.91) . is formed from integer powers 
of L. Moreover, in view of Lemma [Sj^c), the zeroth power LP = In is redundant and can be 
discarded without affecting the algebras of quasiperiodic sets. Thus, for what follows, we define 
the matrix 


C := 


L-2 

L-i 

L 

£2 


c- 

C+ 


G R 


ooxn 


(4.10) 


which is obtained by “stacking” nonzero integer powers of L one underneath the other and is 
partitioned into the submatrices 


Also, we denote by 





'l' 

C~ := 

L-2 

L-i 

£+:= 



Q := Qoo(£), 

Q- := Qoo(£'), Q+ := Qoo(£+) 


(4.11) 


(4.12) 

(4.13) 


the algebras of quasiperiodic subsets of the lattice generated by the matrices (14.101) and 
(14.111) . so that Q~ and Q+ are subalgebras of Q. 


Definition 12 For a given quantized linear {R, L)-system, the algebra Q of C-quasiperiodic 
sets in (14.721) . generated by the matrix C in ( 14.701) . is called an associated algebra. The algebras 
Q~ and of Cr- and -quasiperiodic sets in (14.731) . generated by the matrices Cr and 

in (14.771) . are called backward and forward algebras, respectively. 


In order to clarify the structure of the associated, backward and forward algebras Q, Q and 
, we define, for any a,b eN, the matrix 


[L-n 


c 


a,b • 


L-^ 

L 


Ha 

.Cj 


G 


jg^(aH-6)rixn 


(4.14) 
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which is partitioned into submatrices and given by 


'L-n- 


■ L ' 


r+ ■— 

y f-N ■ — 


L-i 




(4.15) 


for any iV G N. We will now consider the corresponding algebras of quasiperiodic subsets of 
the lattice TIT'". 

Qa,b ■= Q(a+b)ni^a,b), (4.16) 

Qn ■■= QNniC],), Q+ := QNniC^). (4.17) 

These algebras are monotonically increasing with respect to the corresponding subscripts in the 
sense that 

Qa,b C Qa+l,by Qa,b C Qa,b+ly Qa,b C Qa-\-l^b+ly O-N ^ Qn+1i Q-N ^ O-N+l 

for all a,b, N G M. Furthermore, the algebras Q, Q~ and Q+ in (14.121) and (14.131) are repre¬ 
sentable in terms of (14.161) and (14.171) as 

Q= \J Qa,by Q-=\J Qn^ 2^ = U Qn- 

a,b^l N^l N^l 

Lemma 11 The algebras Qa,b, Q~m Q^j in ( 14.761) and ( 14.771) are transformed by the transi¬ 
tion operator T in ( 14.21) and its set-valued inverse T~^ as follows: 

(a) for any a, 6 G N, 

T ^{Qa,b) c Qa-l,b+l, T{Qafi) C Qa+l,b-l] (4.18) 

(b) for any iV G M, 

T-\Qt,) C T{Q],) c (4.19) 

Proof. We will first prove the assertion (a). Suppose a, 6 G N, and let A be a quasiperidoc set 
A := Q{a+b)niGy Ca,b) ^ Qa,b, where G C is a Jordan measurable set. Then application 

of Lemma [10] yields 


T ^(^) — Q(a+b+l)n 




— Q(^a+b+l)n 




(4.20) 


where the sets G-, G+ C are given by 


G+ := {0} X GT 



R-\0), 


G_:=Gx {0} 



L-^R-\0) 


and inherit Jordan measurability from G and R ^(0). The definition of the matrix Ca,b in (14.141) 
implies that 



T^a—1,6+1 


Ga,l 


L-^ = F_ 


C 


a.+l,6—1 


(4.21) 
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where F+, F_ G {0, i|(«+^+i)"x(“+^+i)"- are permutation matrices 


"0 

0 

In 

0 ■ 


Ian 

0 

0 

o' 

0 

I{a—l)n 

0 

0 

F ■= 

0 

0 

0 

In 

In 

0 

0 

0 

? J. _ . 

0 

0 

I{b—l)n 

0 

0 

0 

0 

^bn_ 


0 

In 

0 

0 


Hence, by using the assertions (c) and (d) of Lemma |7l it now follows from (14.201) and (14.211) 
that 


T (A) — Q{^a+h)n{,Hj^^ Ca-l^h+l)^ T{A) — Q(^a+b)niH-, Ca+l,b-l), 


(4.22) 


where the sets 


\v G : 

u 

G T+G'+ for some u E IF \ 

1 

V 


\u G : 

u 

E F_G_ for some u G Z" 

1 

V 

i 


are Jordan measurable. Therefore, the representations (14.221) imply that T~^{A) G Qa-i,b+i 
and T{A) G Qa+i,b-i, thus establishing the inclusions (14.181) in view of the arbitrariness of the 
set A G Qa,b- We will now prove the assertion (b) of the lemma. Suppose N eN and A G Q^, 
that is, 

2l:=QiVn(G,£+) (4.23) 

for a Jordan measurable set G C Then application of Lemma [TOl leads to 


T ^(A) = Q(N+l)n 




where the set H c is given by 


H :={0}xG + 


In 

r+ 


R-\0). 


(4.24) 


(4.25) 


Since 


C 


N. 


L = in view of (I4.15L then it follows from (14.241) and the Jordan measura¬ 


bility of the set H in (14.251) that 


T-\A) = Q^N+i)n{H, G Q++1. (4.26) 

This representation implies the first of the inclusions in (14.191) due to the arbitrariness of A G 
Qlq. The second of the inclusions (14.191) can be established in a similar fashion by using the 
relations 


T{QNn{G, Cj^)) — (5(iV+l)n(G^ X {0} (0), G Qjv+1 

which hold for any N eN and any Jordan measurable set G C which completes the proof 
of the lemma. ■ 
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Theorem 6 The associated, backward and forward algebras Q, Q~ and Q+ in d4.i2D and 
are transformed by the transition operator T in d4.2D and its set-valued inverse T~^ as 
follows: 

(a) the associated algebra Q is invariant under the maps T and T ^ : 

T-\Q) C Q, T{Q) c Q; 

(b) the forward algebra is invariant under the map T~^: 

T-\Q^) c Q+; 

(c) the backward algebra Q is invariant under T: 

T{Q-) C Q-; 

(d) for any A G Q, there exists G N such that T~^{A) G and T^{A) G Q~ for all 
k^N. 

Proof. The assertions (a), (b) and (e) of the theorem follow from Lemma [TT] The assertion 
(d) ean be proved by using the inelusions 

c c Q+, c c Q- 

(whieh hold for all a, 6, c G N and follow from Lemma [TT]) and the property that for any A E Q 
there exist a, 6 G N satisfying A G Qa,b- ■ 

By Theorem[6tb), the transition operator T is measurable with respeet to the forward algebra 
Q+. The assertion (d) of the theorem ean be interpreted as an absorbing property of Q+ with 
respeet to the map and the absorbing property of the baekward algebra Q~ with respeet to 
the transition operator T. 

4.3 Frequency preservation on the forward algebra 

We will need the following enhancement of the nonresonance property for the matrix L of the 
quantized linear (R, L)-system. 

Definition 13 A nonsingular matrix L G is said to be iteratively nonresonant if the cor¬ 
responding matrix C G R°°^” associated with L by M.10\i . is nonresonant in the sense of 
Definition 

Since the matrix L is nonsingular, then the iterative nonresonance property is equivalent to the 
rational independence of the rows of the matrix 


N, 


■‘■n 

L 




Dn(A^H-l)xn 
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for every N e N, where is the matrix given by (14.151) . Therefore, L is iteratively nonres¬ 
onant if and only if so is L~^. Also note that iteratively nonresonant matriees L G do 
exist. Moreover, iteratively resonant matriees L (whieh are not iteratively nonresonant) form a 
set of zero n^-dimensional Lebesgue measure. Indeed, the set of sueh matriees L G ean 
be represented as a eountable union 


U U 


of the following sets 


N 


^N,u,v ■■= \ L^ : det L ^ 0 and ^(L^) V = v\, 


k=l 


where u := 



is an A^n-dimensional veetor partitioned into n-dimensional subveetors 


UN. 

Ui,.. .,un- Eaeh of the sets has zero -dimensional Lebesgue measure whieh ean 

be verified as follows. By assuming, without loss of generality, that \un\ = 1, the Hilbert 
spaee (endowed with the Frobenius inner produet of matriees @) ean be split into the 
orthogonal sum = span(Z) © Z-^ of the one-dimensional subspace spanned by a nonzero 
idempotent matrix Z := unuJj and the corresponding orthogonal complement Z-^ := {Y G 
. T^YiY'^Z) = uJfYuN = 0} wMch is an — 1)-dimensional hyperplane in 
Now, by considering an R"^-valued polynomial h{X) := J2k=i — v for X := Y + XZ, 
with Y G Z-^, it follows that h(Y + XZ) is a polynomial of degree N with respect to A G R, 
with the leading coefficient Z^un = unuJjUn = kat 7 ^ 0 in view of the idempotence of Z. 
Therefore, the following integral with respect to the — 1) -dimensional Lebesgue measure 
over the hyperplane Z-*- vanishes. 


mes„ 2 T 7 v,n,D ^ mes„ 2 |X G R"'^"' : h{X) = O} 

= [ mesi{A G R : hiY + XZ) = 0}dL = 0 
Jz^ 

because the set, involved in the integrand, is finite (consisting of at most N values of A) and, 
thus, has zero one-dimensional Lebesgue measure for every matrix Y G Z^. 

Theorem 7 Suppose the matrix L G of the quantized linear {R, L)-system is iteratively 
nonresonant. Then: 


(a) the associated algebra Q in ( 14.721) consists of frequency measurable subsets of the lattice 

(b) the backward and forward algebras Q and in (14.731) are independent in the sense 
that 

F = F{A)F{B) 

for all A G Q~, B G Q^; 

(c) the transition operator T in ( 14.21) preserves the frequency F on the forward algebra : 

F(r-^(A)) = F{A) for all A G Q+. (4.27) 
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Proof. The assertion (a) of the theorem follows from Theorem lU^a). The assertion (b) is a 
corollary from Theorem H^b) since the matrices C and £+, defined by ( 14 . 111 ) . are nonoverlap¬ 
ping submatrices of the matrix £ given by ( 14 . 101 ) . We will now prove the assertion (c). Suppose 
A E Q'^ is a fixed but otherwise arbitrary set from the forward algebra, and hence, A is repre¬ 
sentable by ( 14 . 231 ) for G M and a Jordan measurable set (7 C [0,1)^". Then, as was obtained 
in the proof of Lemma fTTT bl. the set T~^{A) is given by ( 14 . 261 ) . where H C is a Jordan 

measurable set defined by ( 14 . 251 ) . Therefore, application of Theorem [ 3 ] (under the assumption 
that L is iteratively nonresonant) yields 


F(A) = mes NnG, 

F{T-\A)) = mes^N+i)n , 


(4.28) 

(4.29) 


where V is an arbitrary Lebesgue measurable cell in It will be convenient to use the 

set 

(/;= i?-i(0) X [0,1)^” (4.30) 

which is a cell in in view of Lemma |5];d). From (14.251) and (14.301) . it follows that the 

set on the right-hand side of (14.291) is representable as 


where 


Hence, 


K : = 


V p|(i/ + 


u 

V 


W^:= [0,l)^"fl(£+ 


meS(Ar+i)n 



: uE R ^(0), V E IF^I , 


u + G + Z^^). 

(4.31) 

mes Nn^udu. 

(4.32) 


R-i(O) 


Here, the integrand is identically constant since (14.311) and the above assumption that G C 
[0,1)^” imply that mesMnWu = vnesMnG for all u E M". In combination with the equality 
meSni?“H0) = 1, this reduces the integral in (14.321) to mes(Ar+i)nA' = mesivnt^^, and hence, 
(14.291) takes the form 


F(T-i(2l)) =mes^„G. 


(4.33) 


In view of arbitrariness of the set A E Q'^, comparison of (14.281) with (14.331) establishes (14.271) . 
thus completing the proof of the theorem. ■ 

Theorems [6] and Hare illustrated by Fig.[2l In particular, in view of Theorem |7];b) (under the 
assumption that the matrix L is iteratively nonresonant), F(74) E {0,1} for all A E Qt f) 2^- 
Indeed, any set A E Q~ f) is self-independent in the sense that F(74) = F (^4 f) A) = 
(F(A))^ and hence, F(A) equals either zero or one. This is a version of Kolmogorov’s zero-one 
law [fftll . 

Theorem |^b) and Theorem |3c) show that, under the iterative nonresonance assumption 
on L, the quadruple (Z"^, Q+, F, T) can be regarded as a dynamical system with an invariant 
finitely additive probability measure F, which can be studied from the viewpoint of ergodic 
theory. In particular. Theorem [TT] in the next section establishes a mixing property for this 
quadruple. 

An ergodic theoretic result is provided by Theorem[^below. Note that the transition operator 
T preserves the frequency F on the forward algebra Q+, and this property does not necessarily 
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frequency 
preservation 
and mixing 

Figure 2: The associated algebra Q is invariant under the transition operator T and its set¬ 
valued inverse T~^. The backward algebra Q~ and the forward algebra Q+ are invariant under 
T and T~^, respectively. Under the assumption that the matrix L is iteratively nonresonant, the 
algebra Q is contained by the class S of frequency measurable subsets of The map T is 
frequency preserving and mixing on Q+. The algebra Q" P| Q+ consists of trivial sets whose 
frequencies are either zero or one. 

hold for the associated algebra Q. It turns out that the restriction F|q+ can be extended to a 
finitely additive probability measure on the whole associated algebra Q in such a way that the 
extended measure is preserved under the transition operator. 

Theorem 8 Suppose the matrix L G of the quantized linear {R, L)-system is iteratively 
nonresonant. Then: 

(a) for any element A e Q of the associated algebra Q in ^4.721) . there exists a limit 

F(2l) = lim F(T-^(2l)); (4.34) 

/c^+oo 


(b) the functional F : Q —)■ [0,1] is a finitely additive probability measure on (Z”, Q) and 
satisfies 

¥{T-\A)) = F(2l) ^ F(T(2l)) for all A G Q; (4.35) 

(c) the measures F and F coincide on the forward algebra d4.73D . 

F(2l) = F(y4) for all A G Q+. (4.36) 

Proof. In order to prove the assertion (a) of the theorem, we fix an arbitrary A E Q. The 
absorbing property of the forward algebra Q'^ with respect to T~^ from Theorem [^d) implies 
that there exists N E N such that T~^{A) G for all k N. Hence, from the frequency 
preservation property of the transition operator T on Q'*’, established in Theorem|7]^c), it follows 
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that F(T ^{A}) = F(T ^(A)) for all k ^ N, which implies the convergence (14.341) . In order 
to prove the assertion (b), we note that the property that F is a finitely additive probability 
measure on Q) follows from that of the functional F o T~^ ; Q —)■ [0,1] for any fc G N. 
The equality in (14.351) follows from the definition of F. The inequality in (14.351) can be obtained 
by using the inclusion A C T~^(T{A)) for any A C lA. Indeed, by combining this inclusion 
with the properties of F established above, it follows that 

F{A) ^ F{T-\T{A))) = F{T{A)) 

for any A G Q. Finally, the assertion (c) is proved by noting that (14.361) follows from the 
measurability and frequency preservation properties of the transition operator T with respect to 
the forward algebra Q+ (see Theorem |^b) and Theorem l^c)), which completes the proof of 
the theorem. ■ 


4.4 Independence and uniform distribution of quantization errors 


We will now apply the frequency-based analysis on quasiperiodic subsets of the lattice to 
the deviation of trajectories of the quantized linear {R, L)-system from those of the original 
linear system with a nonsingular matrix L G In one step of the system dynamics, such 

deviation is described by a map E : Z” —)■ i? ^ (0) defined by 

E{x) := Lx — T{x) (4.37) 


More generally, for any k E N, the deviation of trajectories of the systems in k steps of their 
evolution can be expressed as 


N 


L^x-T^{x) = 


k[x 


(4.38) 


k=l 


in terms of maps E^ ■. lA ^ R ^ (0) defined by 

Ek:=EoT^-\ (4.39) 

Definition 14 The map : Z” —)■ i?“^(0), defined by ( 14.371) and (14.391) . is called the kth 
quantization error. 

Also, for any G N, we define a map 8^ '■ "AA ^ (i?“^(0))^ whieh is formed from the 
first N quantization errors as 


8n 


E 


N 


(4.40) 


Lemma 12 For any N e'N, there exists a map qn ■ -E- {R ^(0))^ such that 

(a) the map 8n in ( 14.491) is representable as 

8n = dN ° 8'^, 


(4.41) 


where the matrix C'jif G 


pNnxn 


is given by ( 14.751) .- 


(b) the map is unit periodic with respect to its Nn variables; 


(c) qn is mesNn-continuous; 

(d) is a me.sjqri-preserving bijection of the set {R ^(0))^ onto itself. 
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Proof. We will carry out the proof by induction on N. By recalling (14.21) . it follows that the 
first quantization error Ei = E in (14.371) is representable as 


£i = E = g o L, 

where the map g -.W^ ^ is given by 


(4.42) 


g{u) = u — R{u). (4.43) 

In view of the commutation property (14.11) and Jordan measurability of the set i?“^(0), the map 
g is unit periodic with respect to its n variables and is mesn-continuous. Furthermore, g maps 
the set i?“^(0) identically onto itself and hence, preserves the n-dimensional Lebesgue measure 
on this set. Therefore, the map g in (14.431) indeed satisfies the conditions (a)-(d) of the lemma 
for iV = 1. Now, assume that the assertion of the lemma holds for some G N. Consider the 
next map £n+i ■ ^ in (ICTl) : 


From (14.381) . it follows that 




8n 

En+1 


(4.44) 


N 


T^{x) = L^x-^L^-^E, 




fe=i 


which, in combination with (14.421) and (14.431) . leads to 

En+,{x) := E{T^{x)) = g{L^^^x - E^8 n{x)), 

where the matrix F/v £ is given by 

E^ := [L^ ... L]. 

In view of (14.411) . the representation (14.451) implies that 

En+i = h]\f o 

where the map R~^{0) is defined by 



( 

/ 

Ml 

\ 

\ 


Ml 

hN{u) := g 

Un+ 1 — EjsigN 





, u := 



V 

V 

Ujsf 




Un+1 


(4.45) 


(4.46) 


(4.47) 


where the vector M G ^(iV+l)n 

is partitioned into subvectors ui,... ,un+i G M". By substituting 
(14.461) into (14.441) and recalling (I4.41L it follows that 8n+i = Qn+i ° where the map 

gN+i ■ —)■ (7?“^(0))^+^ is given by 


gN+i{u) ■ — 



(4.48) 
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Note that the map Qn+i is the skew product IfT^fTSlI of the maps and From the properties 
of the map g in (14.431) . it follows that the map in (14.471) is unit periodic with respect to 
its (A^ + l)n variables (which are the entries of the vectors ui,... ,un+i) and is mes(Ar+i)„- 
continuous. Moreover, for any fixed subvectors the map 


h 


N 


\l 


Ml 


Un 


\ 




which is obtained in (14.471) from g by translating the argument, is unit periodic in its n variables, 
mes„-continuous and bijectively maps the set onto itself, preserving the n-dimensional 

Lebesgue measure on this set. Therefore, the map in (14.481) is unit periodic in its (A^ + l)n 
variables, mes(Ar+i)n-continuous and bijectively maps the set onto itself. The 

property that preserves the (A^ + l)n-dimensional Lebesgue measure is established by a 
similar reasoning as for the skew products of measure preserving automorphisms. This com¬ 
pletes the induction step and the proof of the lemma. ■ 


Theorem 9 Suppose the matrix L is iteratively nonresonant. Then the quantization errors 
in ( 14.391) are mutually independent, uniformly distributed on the set R~^ (0) and are measurable 
with respect to the forward algebra Q^. More precisely, for any A^ G N, the map £n in (14.401) 
is mesNn-distributed over the set (i?“^(0))^, and the algebra Ssj.^, generated by this map, 
coincides with the algebra given by ( 14.771) . 

Proof. For a fixed but otherwise arbitrary N G M, let ^fAr : —)■ (7?“^(0))^ be the map 

satisfying the conditions (a)-(d) of Lemma[T^and constructed in the proof of the lemma. Then 
(14.411) implies that 

£]^\B) = {x G Z" : r+x G gf\B)} (4.49) 

for any set B C (77“^(0))^. For what follows, B is assumed to be Jordan measurable. Now, 
since 7?“^(0) is a Jordan measurable cell in M", then, in view of Lemma[5td), the set (7?“^(0))^ 
is a Jordan measurable cell in In combination with the condition (b) of Lemma [T^ this 
implies that gf^{B) is a translation invariant subset of R^"^ which, in view of Lemma [^b), is 
representable as 

gf\B) = G + Z^f (4.50) 

with 

G:=(7?-i(0))^fl<?)^^(i7). (4.51) 

Due to the condition (c) of Lemma [T^ the set G is Jordan measurable. Therefore, from (14.491) 
and (14.501) . it follows that 

£f\B) = Q^r^{G,Ct,), (4.52) 

and hence, this set is £^-quasiperiodic. By Theorem [3] under the iterative nonresonance as¬ 
sumption on L, the set (B) is frequency measurable, with 

F(£:j^^(5)) = mesNnG. (4.53) 
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By recalling (14.511) and the mesTVn-prcserving property of the map (see the condition (d) of 
Lemma [T^. it follows that mesTVnG = mesATni? and hence, (14.531) takes the form 


F(£:^^(S)) = mesNnB. 

Since the last equality holds for an arbitrary Jordan measurable set B C (i?“^(0))^, then the 
map Sn is indeed mesTVn-distributed (that is, uniformly distributed) over the set (i?“^(0))^). 
Finally, by using the property that bijectively maps the set (i?“^(0))^ onto itself (see the 
condition (d) of LemmafT^. it follows that for any Jordan measurable set G C (0))^, there 
exists a Jordan measurable set B C (i?“^(0))^ which represents G by (14.511) . In combination 
with (14.521) . this implies that the inclusion C Q% is, in fact, an equality Ssj^ = which 
completes the proof of the theorem. ■ 


Theorem 10 Suppose the matrix L G of the quantized linear {R, L)-system is iteratively 
nonresonant. Then the algebras Se^, generated by the quantization errors Sk in (14.591) . are 
representable as 

Se, = C Qt for all keN, (4.54) 


and are mutually independent in the sense that for any N and any sets Ak & Se,., ^ ^ k ^ 


N, 



(4.55) 


Proof. The equality in (14.541) is obtained by using the definition (14.391) and the relation Se = 
Qf which was established in the proof of Theorem |9l The inclusion in (14.541) follows from 
Lemma [TTT bf. We will now prove the mutual independence of the algebras generated by the 
quantization errors. To this end, suppose iV G N, and let G Se^. be arbitrary elements of the 
corresponding algebras, representable as 

Afc = Ef\Gk) = Sf\{R-\0))^-^ X Gk), k = l,...,N, (4.56) 

where Gk are Jordan measurable subsets of i?“^(0), and use is made of the maps £k defined by 
(14.401) . Then application of Theorem |9] leads to 

F(Afc) = meSkn{{R~^{0)f~^ x Gk) = meSnGk, (4.57) 

where the rightmost equality also follows from the frequency preservation property of the tran¬ 
sition operator T on the forward algebra Q'^. On the other hand, the intersection of the sets Ak 
in (14.561) is representable as 


N 

[]Ak = Sf\GiX ...xG^) 

k=l 

and hence, by Theorem |9l its frequency takes the form 


F 



N 

mesNn{Gi X ... X Gn) = ]^mesn(3'fc. 

k=l 


(4.58) 
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A comparison of (14.581) with (14.571) leads to (14.551) . which completes the proof of the theorem. 


In addition to extending the results of ifTSl on the asymptotic distribution of roundoff er¬ 
rors, Theorems l9] and [T^ clarify the role of the forward algebra Q'^. More precisely, under the 
assumption that L is iteratively nonresonant, Q'^ is the minimal algebra containing all the mu¬ 
tually independent algebras of frequency measurable quasiperiodic subsets of the lattice, 
generated by the quantization errors for k E N. Equivalently, the forward algebra describes 
events which pertain to the deviation of positive semitrajectories of the quantized linear sys¬ 
tem from those of the original system over finite time intervals (which will be considered in 
Section |431) . 

The following theorem is a corollary from Theorem |9] and develops an ergodic theoretic 
point of view for the quantized linear (R, L)-system, for which, as mentioned before, the re¬ 
striction of the frequency F to the forward algebra Q~^ is an invariant finitely additive probability 
measure. 

Theorem 11 Suppose the matrix L E is iteratively nonresonant. Then the quadruple 

[IT., Q'*', F, T) satisfies the mixing property 

lim FfT-^(A)n5') =F(A)F(5), A,BeQ^. (4.59) 


Proof. Let A, B E be fixed but otherwise arbitrary elements of the forward algebra. Then 
there exists N eN such that A,Be where the algebra is defined by (14.171) . Hence, in 
view of Theorem m the sets A and B are representable as 

A = S^\G), B = S],\H) (4.60) 

for some Jordan measurable sets G, H C . From the definitions (14.391) and (14.401) . it 

follows that 

Ek+l 

£noT'^= : 

Ek+N 

for any A; ^ 0. In combination with the first of the equalities in (14.601) . this implies that the set 


T-\A) = {Sj^oT^)-\G) 


|x e Z” : 


Ek+fix) 

Ek+N{x) 


E G 


belongs to the algebra generated by the quantization errors Ek+i: ■ ■ ■, Ek+N- Since 

Ek+i 


for every k ^ N, the maps £n and 


are formed from nonoverlapping sets of quantiza¬ 


tion errors, they induce mutually inc 
from (14.601) . it follows that 


Ek+N 

ependent algebras. Hence, by recalling the second equality 


F (^T-^(A)p|5) = F(T-'=(A))F(5) = F(A)F(5). 
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Since the last two equalities hold for any k ^ N, they imply the eonvergenee (14.591) . thus 
eompleting the proof. ■ 

Note that the proof of Theorem [H] establishes a stronger result: under the assumption that 
L is iteratively nonresonant, the algebras and are mutually independent for all 

N E N and k ^ N, and hence, by the version of Kolmogorov’s zero-one law of Section |431 the 
algebra Q'^ consists of trivial subsets of the lattice whose frequencies are either 

zero or one. 


4.5 Distribution of the deviation of positive semitrajectories of the quan¬ 
tized linear and supporting systems 

For what follows, we define an affine map T* : M” —)■ M" by 

T^{x) := Lx — fi, (4.61) 

where 

/i := / udu (4.62) 

is the mean veetor of the uniform probability distribution over the cell i?“^(0). 

Definition 15 For the quantized linear [R, L)-system, the dynamical system in with the 
transition operator T* given by ( 14.(571) is referred to as the supporting system. 

For every G N, we introduee a deviation map 5^ '■ Z" —)■ M", whieh deseribes the 
deviation between the A^th iterates of the quantized and supporting systems as 


N 


k=l 


N 


5^{x) := T^{x) - T^{x) = - d)- 

Also, we define another map —)■ ML by 

(„(x) ~x- T.-'^(T'^(x)) = L-'=(Et(x) - IX). 

k=l 

The maps 5^ and are related to eaeh other by a linear transformation 

6n = 

and satisfy reeurrence equations 

5m = L5m-i + Em — p, = ^n-i + L ^{Em — p) 


(4.63) 


(4.64) 


(4.65) 


(4.66) 


for all 24 G N, with initial conditions = 0. 

For every A^ G M, the map ^m in (14.641) is an affine funetion of the first N quantization 
errors. Henee, in view of Lemma|4]and Theorem|9]for any iteratively nonresonant matrix L, a 
map Sat : —>■ defined by 


xN .— 


I 

Uji * 


1 

■ ^ ■■ 



L-^ 

L-i 


0 

L-n 



1 

1 . 


1 
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is uniformly distributed over the set 


L-i 

L-i 


0 

L-n 




with constant probability density function (PDF) | det L\ 2 with respect to mesTVn- More¬ 
over, the corresponding algebra coincides with the algebra Q'^ defined by (14.171) . 

Therefore, the maps 6^ and can be regarded as M’^-valued random vectors on the prob¬ 
ability space (Z”, Q+, F). From (14.661) and the mutual independence and uniform distribution 
of the quantization errors, it follows that the sequence {Sn)n^i has the structure of a homo¬ 
geneous Markov chain, whilst {^n)n^i is a non-homogeneous Markov chain with independent 
increments (whose average values are zero since A{Ek) = /r in view of (14.621) 1. We will use 
these properties in Section 15.21 in order to establish a functional central limit theorem for the 
sequence for a class of quantized linear systems. 


4.6 Distribution of the deviation of negative semitrajectories of the quan¬ 
tized linear and supporting systems 

For any fc G N, the preimage 

T-\x) := {yez^-. T\y) = x} 

of a point x G Z" is a finite (possibly empty) subset of the lattice Z". The set-valued sequence 
(T-"(x))fe^i describes a negative semitrajectory of the quantized linear (i?, L)-system, and its 
fragment 

T-\x),...,T-^{x) (4.67) 

is a basin of attraction of depth N for the point x. The aim of this subsection is to study 
the deviation of negative semitrajectories of the quantized linear (72, L)-system from those of 
the supporting system introduced in the previous section. To this end, we define yet another 
dynamical system on IP with the transition operator T -.IP ^ IP given by 

f{x) :=72(/i + r-'(x)), (4.68) 

where y is the mean vector from (14.621) . and T* is the transition operator of the supporting 
system in (14.611) with the inverse 

T“^(x) = L~^{x + y). 

Therefore, the map T in (14.681) is representable as 

f = RoL-\ (4.69) 

where 72 : M" —)• Z" is a quantizer given by 

72(m) := R{u -f {In + L~^)y)- (4.70) 

Definition 16 For the quantized linear {R, L)-system, the quantized linear {R, L~^)-system 
with the quantizer 72 in ( 14.701) . is called the compensating system. 
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Note that the map {R, L) i—)■ (i?, L~^) is an involution in the sense that the eompensating 
system for the quantized linear (i?, L“^)-system eoineides with the quantized linear (i?, L)- 
system. ^ ^ 

We will now define quantization errors : IT- ^ for the eompensating system 

in a similar fashion to the quantization errors for the quantized linear (i?, L)-system in 
Seetion l4^ ^ _ _ 

Ek:=EoT^-\ = 1,2,3,..., (4.71) 

where 

E{x) := L~^x — T{x). (4.72) 

The quantization errors Ek of the eompensating system take values in the eell i?“^(0) whieh, in 
view of (14.701) . is given by 

^-1(0) =i?-i(0)-(4 + L-i)/i. (4.73) 

Theorem 12 Suppose the matrix L G is iteratively nonresonant. Then: 

(a) the quantization errors 77^ of the compensating system in ( 14.771) are mutually independent, 
uniformly distributed on the cell 77“^ (0) in (14.731) and are measurable with respect to 
the backward algebra Q~ given by (14.731) . More precisely, for any N E N, the map 
^ defined by 


£n '■= 



(4.74) 


is mesNn-distributed, and the algebra generated by this map, coincides with the 
algebra Qf defined in (14.771) .- 

(b) the algebras, generated by the quantization errors of the compensating system, are rep¬ 
resentable as C Of for all k gN, and are mutually independent. 


Proof. It was remarked in Seetion 14.31 that if the matrix L is iteratively nonresonant then so 
is its inverse L~^. Therefore, the assertions of the theorem can be established by applying 
Theorems l9l and [TOl to the compensating system. ■ 

Now, for any G N, we define a set-valued map : IT fF' by using the transition 
operator T of the compensating system as 

S^(a:) := T-^{x) - f^{x). (4.75) 

This allows the preimage T' ^ {x) of a point a: G to be represented as the translation of the 
set Sjv(a;) C Id by the vector T^{x) G Id: 

T-^{x) = f^{x) + (4.76) 

In particular, for = 0, the definition (14.751) implies that the set So(x) is a singleton consisting 
of the origin: 

So(a:) = {0}. (4.77) 
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In view of this initial condition, it follows by induction that the sets in (14.751) satisfy the 

recurrence equation 


S^(a;) =Z’^p|(Ejv(a;)+L-^(S^_i(a;)+i?-'(0))), NeN. (4.78) 

Indeed, by using (14.751) . (14.761) and recalling the first quantization error E from (14.721) . it follows 
that 

= (z- f| (^L-^f^-\x) - f^{x) + l-\j:n_,{x) + i(-'(0)))) 

= (z^ fl (E{f^-\x)) + L-\Zn_i{x) + i?-'(0)))) 

= Z" f| + L-\Zn_,{x) + R-\0))) 

for any iV G M. For what follows, we add the empty subset of the lattice to the class 11 in (12.11) 
and denote by 

n := {21 C Z" : #2l < +cx)} = U 1J{0} (4.79) 

the class of finite subsets of Z”. Also, we define a function A : 11^ —)■ [0,1] which maps a pair 
of such sets A and B to 

A(B I A) = mes„ |m G R-^(O) : Z" p| (m + L-^(A + 7C^(0))) = . (4.80) 

Note that Xlseff I ^) = 1 for any A G 11, and hence, the function A(- | A) ; 11 ^ [0,1] 
describes a conditional probability mass function on the class 11 (which is a denumerable set). 

Theorem 13 Suppose the matrix L is iteratively nonresonant. Then the sequence of the set¬ 
valued maps : Z"^ —)■ n defined by ( 14.751) is a homogeneous Markov chain with respect to 
the filtration formed from the algebras in (14. 1 71) with the transition kernel A(- | ■) in 

( l4.igOD . More precisely, for any A G N and any finite subsets Bi,of the lattice Z” (that 
is, elements of the class 11 in ( 14.791) ]. 

N 

fl Z-^\Bk) = {x G Z" : Sfc(a;) = 5, for all fc = 1,..., A} G 

k=l 

and 

( N \ ^ 

k=l J k=l 

with Bq := {0}. 
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Proof. From (14.771) and (14.781) . it follows by induction that for any iV G N, the map in 
(14.751) is representable as 


Sat — Pn ° S-N- 


(4.81) 


Here, 8n is the map defined by (14.741) and pN '■ {R —)■ IT is a mesAr^-eontinuous map 

satisfying the reeurrence equation 


p^(u) = Z"f| 


( ( 

Ml 

\ 

\ 

\ 

un + ^ Pn-1 



+ R-\0) 


\ \ 

Un-1 



/ 


(4.82) 


where the vector u 



E (R ^(0))^ is partitioned into n-dimensional subveetors 


lUN. 

Ml,... ,MAr G i?“^(0), and the initial eondition pQ = {0} is used. The assertions of the the¬ 
orem ean now be obtained from the representation (14.811) and the reeurrenee relation (14.821) by 
using Lemma |4] and Theorem [T^ ■ 


The relation (14.761) implies that 


T-^ix) - T8^{x) = Siv(a;) - 5iv(a;), 


(4.83) 


where the maps 6^ ■. IR ^ are defined by 


N 

Sn := + L-V) 

k=l 

and satisfy the reeurrenee equation 

5]\f = L + E]\f + L ^/i, N = 1,2,3,..., 

with the initial condition (5o = 0. This equation is similar to the first of the reeurrence relations 
(14361) . 

Therefore, the importanee of Theorems [I^ and [l3] eonsists in representing the deviation 
(14.831) as the Minkowski sum of two sequenees driven by mutually independent and uniformly 
distributed quantization errors of the eompensating system: the n-dimensional sequenee 5n{x) 
and the set-valued sequenee Sat (x) whieh are homogeneous Markov ehains with respeet to the 
filtration formed from the algebras (14.171) . 

This elarifies the role of the backward algebra Q~ defined by (14.131) : while the forward 
algebra Q+ is formed from events whieh pertain to the deviation of positive semitrajeetories 
of the quantized linear (R, L)-system and the supporting system (see Seetion 1431) . the back¬ 
ward algebra Q~ deseribes events related to the deviation of negative semitrajectories of these 
systems. 

Now, for any /c G N, we define a funetion Uk : 1R ^ 1^+ by 

z/fc(x) := #r-'^(a;) = #Sfc(a;) (4.84) 


The eorresponding map := 


1^1 


un 


the basin of attraetion (14.671) . From T 


: lA -E- describes the cardinality structure of 


leorem [13] (under the assumption that the matrix L is 
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iteratively nonresonant), it follows that for any N G N, the map A/iv is distributed and C 
Q~n- 

Theorem 14 Suppose the matrix L is iteratively nonresonant. Then the following functions, 
obtained by rescaling H.84^ as 

un ■= I det L\^un : Z” —>■ M+, (4.85) 

form a martingale ([74|/ with respect to the filtration {Qf}N^i of the algebras M.lTil in the 
sense that 

A(i/jv+i2A) = for all A e N, A e Q~f. (4.86) 


Proof. Sinee, by Theorem[T3l the sequenee (S^r) jv^i of the set-valued maps (14.751) is a homo¬ 
geneous Markov ehain with respeet to the filtration {Q'f}N^i with the transition probabilities 
(14.801) . the statement of Theorem [T4l will be proved if we show that 

V A(5 I A)fiB = , , for all AeU. (4.87) 

detL 

Ben 


To this end, we will first prove that if V, (7 C M” are Lebesgue measurable and 1/ is a eell, then 

[ #(^Z"p|(M + G))dM = mes^C;. (4.88) 


'V 


Indeed, 


[ #(z^n{u + G))du= [ y^I^+G{z)du 
Jv ^ ^ Jv 

= [ Ig{z- u)du = ^ mes„(G' p|( 2 ; - V)) 

= meSn^f^Pl (Z"" — fo) j = raeSnG. 

A eombination of (14.801) and (14.881) allows the sum on the left-hand side of (14.871) to be repre¬ 
sented as 


Va(BIA)#B= [ fi(z-n{u + L-\A + R-\t))))yu 

mes„(A i?“^(0)) 


Ben 


= mes„(L {A + R (0))) = 


det LI 


(4.89) 


Since i?“^(0) is a cell in M" and A is a finite subset of Z"^, then mes„(A -f i?“^(0)) = fiA. 
Hence, the right-hand side of (14.891) leads to that of (14.871) . which completes the proof of the 
theorem. ■ 


The martingale property (14.861) of the rescaled functions (14.851) established by Theorem [T4l 
implies that (under the assumption that L is iteratively nonresonant) the average values of the 
functions (14.841) can be computed as A(z/fc) = | det L|“^ for all k &N. 

In order to provide another corollary from Theorem[T3l we will now consider the set {TA) 
of those points of Z” which are reachable in N iterates of the transition operator T. These 
reachability sets satisfy 

yiv+i(^n) ^ = {xez^ ■. un{x) > 0} for all N en. (4.90) 
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Theorem 15 Suppose the matrix L is iteratively nonresonant. Then for any N ^ N, the fre¬ 
quency of the Nth reachability set in ^4,901) can be computed as 


FiT^iZ^)) = H^m) 

in terms of functions i^Ar : 11 —)■ [0,1] which are governed by the following linear recurrence 
equation 

Sen 

for all A E A, with the initial condition Hq = 1, where A(- | ■) is the Markov transition kernel 
from \4. ^01) . 

Proof. The assertion of the theorem follows from Theorem [13] and from the property that 0 
is an absorbing state of the set-valued Markov ehain (SAr)jv^i in its state spaee If defined by 
dHH). " ■ 

From Theorem [T5l it follows that, under the assumption that the matrix L is iteratively 
nonresonant, 

F(T(Z")) = Hf{0}) = 1 - A(0 I {0}), (4.91) 

where ^ 

A(0 I {0}) = mes^jw G R-\0) : Z" p| (m + L-^R-\0)) = 0| (4.92) 

is the frequeney of the set Z^ \ T(Z") of “holes” in Z"- (not reaehable for T) whieh quantifies 
nonsurjeetivity of the transition operator T. Moreover, sinee the relation (14.911) involves only the 
first iterate of T, it remains valid if the nonsingular matrix L is resonant (that is, not neeessarily 
iteratively nonresonant), with T(Z") being a frequeney measurable L“Aquasiperiodie subset 
of Z”. 


5 Quantized linear systems: neutral case 

5.1 The class of quantized linear systems being considered 

We will now eonsider a partieular elass of quantized linear (i?, L)-systems on the integer lattiee 
Z'^ of even dimension n := 2r, with the matrix L G being similar to an orthogonal matrix: 


L:=J2UkJkVk. (5.1) 

k=l 


Here, Uk G and 14 ^ are the bloeks of a nonsingular matrix U G and its 
inverse, 

>r 


U :=[Ui ... Ur] , U-^ : = 


(5.2) 


K- 


and 


Jk ■ 


COS 6 k 

sin dk 


— sin dk 
cos 9k 


(5.3) 
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are the matriees of planar rotation by angles 6k G (0,27r). The matrix L, given by (15.11) - 
(15.31) . has r two-dimensional invariant subspaces For what follows, we assume that the 

spectrum of L is nondegenerate, that is, its eigenvalues 1 ^ k ^ r, are all pairwise 

different. 

Note that, for any given nonsingular matrix U G those vectors 6 := {dk)i^k^r ^ 

(0, 27r)'’, for which the corresponding matrix L in (I5.1I) - (I5.3I) is iteratively resonant, form a set 
of r-dimensional Lebesgue measure zero. Indeed, the set of such vectors 6 is representable as a 
countable union 

1^_J 1^_J ^N,u,v 


of the sets 


where 


'<JN,u,v ■— '^N,u,v ^ ^ '^N,u,vy 


N 


- ■{ «» 6 ( 0 , 27r) : Jk'-Ului = Utv 

£=1 


Here, the vector u is partitioned into n-dimensional subvectors mi, ..., uat. For any iV e N, 
u G \ {0} and u G there exists at least one value of the index fc = 1, ..., r such that the 
set ^ consists of real roots of a nonconstant trigonometric polynomial and hence, is finite. 
Therefore, each of the sets Qn,u,v has zero r-dimensional Lebesgue measure and so does their 
countable union. 

Hence, for any given nonsingular matrix U G the matrix L, described by (15.11) - 

(15.31) . is iteratively nonresonant and has a nondegenerate spectrum for mesr-almost all vectors 
6 G (0, 27r)^ of rotation angles. 

Since the eigenvalues of the matrix L have unit modulus (and hence, the linear dynamical 
system x ^ Lxis marginally stable), the corresponding quantized linear (R, L)-system will be 
called neutral. 


5.2 A functional central limit theorem for the deviation of positive semi¬ 
trajectories of the quantized linear and supporting systems 

In what follows, we will need a technical lemma which is given below for completeness of 
exposition. 


Lemma 13 Suppose the matrix L G given by (I5.i|] -( |53|] . has a nondegenerate spectrum. 
Then, for any real positive definite symmetric matrix 'F of order n, the following convergence 
holds: 


$ := lim 

A^^ + CO 





(5.4) 


where 

- fwvRvf). 


(5.5) 
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Proof. For any N let 


N 






(5.6) 


£=1 


denote the matrix whose eonvergenee is considered in (15.41) . By using (I5.1I) - (I5.2I) . this matrix 
can be represented as 

Vfifjivl, (5.7) 

j,k=l 

where 


N 




e=i 


and 


vh,,, := V,^>V^. 

Each of the planar rotation matrices in (15.31) is representable as 


Jk = W 


Qidk 0 

0 


W*, 


(5.8) 

(5.9) 

(5.10) 


where (•)* denotes the complex conjugate transpose, and 


W := 


1 1 


\/2 i 

Substitution of (15.101) into (15.81) leads to 

where 0 denotes the Hadamard product of matrices [[^, and 


N 


^(JV) _ A V 


e=i 







-1 

'Ip 

_ 1 


(5.11) 


(5.12) 


(5.13) 


Under the assumption that the spectrum of the matrix L is nondegenerate, the following con¬ 
vergence holds: 

I 2 for j = k 
0 for j ^ k 

A straightforward verification leads to the property of the matrix W in (15.1 II) that 


lim = 

A^^+co 


(5.14) 


TrF 

W (I 2 O (W*FW)) W* = - 12 


(5.15) 


for any real symmetric (2 x 2)-matrix F. Now, from (15.71) and (I5.12I) - (I5.15I) . it follows that 




(5.16) 


k=l 


In view of (15.51) . (15.61) and (15.91) . the convergence (15.161) is equivalent to (15.41) . which completes 
the proof of the lemma. ■ 
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Note that the matrix <h, given by (15.41) and (15.51) . is symmetrie, positive definite and satisfies 
= <h. The latter property implies the orthogonality of the matrix where 

\/$ G I^nxn ^ necessarily symmetric) matrix square root of $ (satisfying \/$(\/$)"'" = 

$) with the inverse := (\/$)“^. It will be convenient to use the following matrix square 
root of $: 

\/$ := ^ [dif/i ... arllr] (5.17) 

v2 


which satisfies 




Ji 


0 


= J. 


0 Jr 


(5.18) 


Now, let the matrix 4/ G in Lemma [13] be the covariance matrix of the uniform 

distribution over the set i?“^(0): 


'^:= {u — ^){u — fi)^du, (5.19) 

Jr-^{0) 

where /r is the mean vector from (14.621) . The matrix 4/ is symmetric and positive definite. Its 
nonsingularity follows from the property that i?“^(0) is a cell in M" and hence, can not be 
contained by an {n — 1)-dimensional hyperplane. 

By using the matrix $ G expressed through (15.41) in terms of the matrix 4/ in (15.191) . 
and recalling (14.641) . we will now define, for arbitrary iV G N and t G [0,1], a map 
R^hy 

:= - ^Livq)) • (5.20) 

Here, [■] denotes the ceiling function, and t plays the role of a continuous time variable. For 
any fixed iV G N and x G Z"^, the map Wn,,{x) : [0,1] —)■ M” is a continuous piecewise linear 
function which satisfies 

k = 0,...,N. 

Hence, the map in (15.201) generates a map Wn : ^ C which maps a point x G Z” to 

the function (x) belonging to the complete separable metric space 

C':=C'o([0,l],M") (5.21) 


of continuous functions / : [0,1] — )■ M” satisfying /(O) = 0. The space C is endowed in a 
standard fashion with the uniform norm ||/|| := max^gjo,!] \f{t)\ and the corresponding Borel 
a-algebra Be- 

Under the assumptions that the matrix L is iteratively nonresonant and has a nondegenerate 
spectrum, the family {W^^t ■ f G [0,1]} of the maps (15.201) can be regarded, for any given 
G N, as an M"^-valued random process on the probability space {Z^, Q~^, F). This follows 
from the property (which can be obtained by using Theorem |9|) that for any p G N and any reals 
0 < < ... < fp ^ 1, the map 


U 




WN,t 




Z" ^ 


(5.22) 
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is distributed with a countably additive probability measure -Dtv,*!,...,* on (M^”, (with 
the (T-algebra of Borel subsets of W^), and the algebra induced by this map is contained by 
the algebra given by (14.171) . The probability measures DN,ti,...,tp satisfy the compatibility 
conditions of A.N.Kolmogorov lfT4l and are the finite-dimensional distributions of the process 
(15:201) . 

On the other hand, under the above mentioned assumptions on matrix L, for any iV G N, the 
map Wn —)■ C, specified by by (15.201) . is distributed with a countably additive probability 

measure on {C, Be), and the algebra induced by this map is again contained by Q'^. There¬ 
fore, the map can be regarded as a random element on the probability space Q+, F) 
with values in the metric space C. The relationship between the probability measures Dj^ and 
DN,ti,...,tp is described in terms of the cylinder sets as 


DN,t,,...,tp{B)=D^{feC: 


fiti) 

fitp) 


G B 


for all B eBP'^. _ 

We will now establish a functional central limit theorem for the map Wp^, that is, the weak 
convergence of the probability measure Dn sls N —)■ -|-cxd, to the countably additive prob¬ 
ability measure W on {C,Bc) corresponding to the n-dimensional standard Wiener process. 
This will be carried out using the standard scheme [HI: first, the weak convergence of the 
finite-dimensional distributions DN^ti,...,tp of the random process (15.201) to those of the standard 
Wiener process will be proved; then it will be shown that the sequence of probability measures 
(Dpf) jv^i satisfies the Prokhorov criterion of relative compactness in the topology of weak con¬ 
vergence of probability measures on {C, Be)- 


Lemma 14 Suppose the matrix L in ( l5.iD - (l5.5D is iteratively nonresonant and has a nondegen¬ 
erate spectrum. Then the finite-dimensional distributions of the random process (15.201) weakly 
converge to those of the n-dimensional standard Wiener process as N ^ -|-oo. More precisely, 
for any p eN, any reals 

tQ:=0 <ti < ... <tp^l (5.23) 

and any Jordan measurable set B C the distribution of the map UN,ti,...,tp in ( 15.221) satisfies 

p 

hm F{x G Z" : UN,tu...,tp{x) G B] ={27i)-p'^B TT(4 _ 4_,)-/2 

N^+oo j J.J. 

k=l 

(5.24) 


where Mq := 0. 


Proof. Let the numbers p E N and ti,... ,tp in (15.231) be fixed but otherwise arbitrary. Then, 
for any N ^ distribution of the map UNfy,...,tp in ( |5.22D is 

absolutely continuous with respect to meSp„. Hence, for any such N, the following equality 
holds for any Jordan measurable set B C 

F{xEZ^: UNM,...,tp{x) EB}= DNfy,...,tp{B). 
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Therefore, the assertion of the lemma will be proved if we establish the weak eonvergenee 
of the probability measure to the joint distribution of the values of the n-dimensional 

standard Wiener proeess at the moments of time (15.231) . In view of the property that the sequenee 
of the maps (14.641) has independent inerements, this task reduces to proving the weak 
convergence of the distribution Dj^f of the map 


N 






-)■ 


(5.25) 


to the n-dimensional Gaussian measure D with zero mean vector and the identity covariance 
matrix. The probability measure can be uniquely recovered from its characteristic function 
—)■ C defined by lfT4l 


(Pn{u) : = 

From Theorem|9]and (15.251) . it follows that 

N 

^n{u )= 

k=l 




(5.26) 


Vn 


($-V2^-fc)T 


U 


(5.27) 


where ^ C is the characteristic function for the uniform distribution over the set 

7?-i(0) -/i: 


: = 




W. 


Ir-Ho) 


(5.28) 


Note that the function (15.281) is infinitely differentiable with respect to n G and its asymptotic 
behaviour near the origin is 


(p{v) = 1 — -||u||| + o(|u|^) as u —)■ 0, 


(5.29) 


where ||u||^ := denotes the Euclidean norm in associated with the covariance 

matrix 4/ in (15.191) . By using (15.291) . (15.271) and Lemma [T3l it follows that 


N 


lim InifNiu) =-- lim 
—2 W—>-+oo \ N ^ ^ 


-k\T^ 


N^-\-oo 


U 


i: 


k=l 


N 




\U\ 


k=l 




holds for any u G M", and hence, lim 7 v^.+oo y’Af(w) = e The latter limit implies that Djq 
converges weakly to the standard normal distribution i9 in M" as iV —)■ +oo, and the lemma is 
proved. ■ 

In order to formulate the functional central limit theorem below, we will now define the 
subalgebra 

Be := js G Be : Wi^dB) = 0 and Dn^OB) = 0 for all sufficiently large N G n| (5.30) 

of the algebra of W-continuous Borel subsets of the metric space (15.211) (as before, W denotes 
the Wiener measure on{C,Be))- 
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Theorem 16 Suppose the matrix L in 0.iD - rt5.3D is iteratively nonresonant and has a nonde¬ 
generate spectrum. Then, for any B G Be, 


lim YixelT : Wnix) E b] = W(B). 

N—>-+oo L J 


(5.31) 


Proof. For any element B G Be of the algebra (15.301) and for all sufficiently large N E N, the 
set W^^{B) = |x G Z" : Wn{x) E i?| is frequency measurable, and its frequency is given 

by F{W^^{B)) = Djsi{B). Since every B E Be is a W-continuous Borel subset of C, the 
relation (15.311) will follow from the weak convergence of the measure DArtoWasiV —)■ +cxo. 
Therefore, in view of Lemma [T4l it remains to show that the sequence satisfies the 

Prokhorov criterion of relative compactness in the topology of weak convergence of probability 
measures on {C,Be)- The latter task reduces to verifying the following condition (see, for 
example, [[B Theorem 8.4]): 


inf f lim sup max F < a; G : 

-^>1 V N^+oo OsSjsJAf ( 


max j (x) — Wj,, k (x) 

j^k^N ^ 



0. (5.32) 


To this end, let iV G N, 0 ^ ^ iV and A > 1 be fixed but otherwise arbitrary. Since, 

under the assumption that L is iteratively nonresonant, the sequence of the maps (14.641) has 
mutually independent increments with zero average value, then a multidimensional version of 
[[B inequality (10.7)] leads to 


F <; a; G Z”' : max 

j^k^N 


= F a: G Z" 


Wj,^±{x) - Wj,f_^{x) 


max 


F 2F X G Z" 


^/N 
1 


^ A 


y/N 


^_i/2 ^ L-^(E,{x) - /i) > A 

e=j+i 
N 

L-\E^{x) -p) ^ a - \/2 !> . 

£=j+l 


(5.33) 


By applying the Chebyshev inequality [[T4ll to the right-hand side of (15.331) . it follows that, for 
any A > \/2, 


F X G Z" 


sfN 


N 






iV2(A- \/2)4 


e=j+i 


A 


^X-y/2 


N 


^_l/2 ^ 

e=j+i 


(5.34) 


The fourth moment on the right-hand side of (15.34!) can be estimated by using the mutual in¬ 
dependence and uniform distribution of the quantization errors over the set i?“^(0). From the 
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orthogonality of the matrix (15.181) and from (15.191) . it follows that 


A 


N 


L-\Ei 

£=j+l 



=A 


N 


i=j+i 


= {N-3) [ 

+ (A-j)(iV-J-l)(TrT)2 


+ 2 Tr 

j+l^£y£m^N 

I fx)fdu +N^K, 


(5.35) 


where the parameter K := (TrT)^ + 2Tr(T^) is assoeiated with an auxiliary matrix T := 
$-i/ 2 ^(<|)-i/ 2 ^T relations (I5.33I) - (I5.35I) lead to the inequality 


limsup max F < x G : max 

N^+oo I j^k^N 


IF. 


(x) - Wj^k{x) 



2K 

(A^Vlr 


(5.36) 


whieh holds for any A > \/2. Since the right-hand side of (15.361) is o(A“^) as A ^ -|-cxo, this 
inequality establishes (15.321) . which, in view of the above discussion, implies the assertion of 
the theorem. ■ 


Theorem 17 Suppose the matrix L in ( l5.iD - (l5.3l) is iteratively nonresonant and has a nonde¬ 
generate spectrum. Then the distribution of the map 

^ M” (5.37) 

VN 

(where use is made of A4.63\l . (15.41) and ( l5.i9D ) weakly converges to the n-dimensional Gaus¬ 
sian measure with the zero mean vector and identity covariance matrix as N —)■ -|-oo. More 
precisely, for any Jordan measurable set B C M”, 

lim fIxgZ": 5n{x) e b \ = {2Tr)-^/^ [ e'^dw. (5.38) 

N^+cxy ( VA J Jb 


Proof. The linear bijection (14.651) between the maps (14.631) and (14.641) allows the map (15.371) 
to be represented as 

^$- 1 / 2 ^^ = ^ jNy^ ( 5 . 39 ) 

VA VA 

where use is made of (15.181) and (15.251) . Therefore, the map in (15.371) is iA^r o J“'^-distributed, 
where is the distribution of the map lFv,i. Hence, for any Jordan measurable B cM."', the 
frequency of the set on the left-hand side of (15.381) coincides with B). Therefore, the 

assertion of the theorem can be proved by establishing the weak convergence of the probability 
measure Djq o J~^ to the ?7,-dimensional Gaussian measure D with the zero mean vector and 
identity covariance matrix as A —)■ -|-cxd. For this purpose, we will use the weak convergence of 
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Dj^ to D as N ^ +00 (which was established in the proof of Lemma [T4l). and the invariance 
of the probability measure D with respeet to orthogonal transformations (sinee standard normal 
PDFs are isotropic). With any given bounded Lipschitz eontinuous funetion / : M” —)■ M, we 
assoeiate the funetions 

/tv := / o NeN. (5.40) 

Due to orthogonality of the matrix J in (15.181) . the functions /at have the same Lipsehitz con¬ 
tinuity constant as the function /. Therefore, in view of the Arzela-Ascoli theorem, for any 
£ > 0, the restrictions of these funetions to a ball ;= {m G : |m| ^ form a totally 
bounded set 

(5.41) 

in the Banaeh space 0(3^, M) of eontinuous funetions on the ball. Henee, there exist G N, 
a surjeetive map : N —)■ {1, 2,..., iVg} and a finite £-net 

{/f : 1 ^ k ^ N,} G C{B,,R) (5.42) 

for the set (15.411) satisfying 

sup |/tv(«) - fu]N){u)\ ^ (5.43) 

AT^l, ueBe 

Now, reealling (15.391) and (15.401) . consider the following integrals 


A / 


y/N 




N - 


/(m)(L>jv o J ){du) 


= / fN{u)DN{du) + 


'Be 


A-Be 


fN{u)DN{du) 


(5.44) 


In view of the weak convergence of to D, the rightmost integral in (15.441) ean be bounded 
asymptotieally as 


lim sup 

W—>-+oo 



fN{u)DN{du) 




(5.45) 


where ||/|| := sup^g^n |/(m)| is the uniform norm of /. A combination of the same weak 
eonvergenee with (15.431) leads to 


liminf / fN{u)DN(du) ^ min / flf\u)D{du) — e. 

N^+oo Jr l^k^Ne Ja 


(5.46) 


On the other hand, from (15.431) and the rotational invariance of the Gaussian measure D, it 
follows that 


mm 

l<k^Np 


fjf\u)D{du) ^ min f fN{u)D{du) — e = f f{u)D{du)—e. (5.47) 

JBe JBe 

A combination of (l5]4^ with (l5]47]) yields lim inf AT^+oo/g fN{u)DM{du)^ /g f{u)D{du) — 
2e. Henee, in view of (15.451) . it follows that 

liminf/' fiq{u)Diq{du)^ f f{u)D{du)-2(e + {l-D{B^))\\f\\). 

N^+oo V / 
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The arbitrariness of e > 0 in the latter inequality and the property lirue^+o D{B^) = 1 imply 
that 

liminf f f]\f{u)D]\[{du) ^ f f{u)D{du). 

N^ + OO J-^ri J-^n 

A similar reasoning leads to the inequality 


lim sup 

A^^ + OO 


fN{u)DN{du) ^ 


f{u)D{du). 


A combination of the last two inequalities leads to the following convergence for the left-hand 
side of (15.441) : 


lim 


N^-\-oo 


f{u){DNoJ ^){du) = [ f{u)D{du). 

7r" 


(5.48) 


Since / : R” —)■ M is an arbitrary bounded Lipschitz continuous function, then, in view of the 
well-known criterion [fT4ll of the weak convergence of probability measures, (15.481) implies that 
Djsf o J~^ weakly converges to D as N ^ -l-oo, whence the assertion of the theorem follows. 


Theorem 18 Suppose the matrix L in ( I5.il) - (I5.5|] is iteratively nonresonant and has a nonde¬ 
generate spectrum. Then for any Jordan measurable set B C Rl, 


lim F < X G 

A^^ + CO 


[Y Jmaxi^k^N \Vj 6 k{x)\ 

- y n[ 

/ x{uk)dui X ... X dUr, 
k=l 


G B 




(5.49) 


where use is made of M.63\) . ( 15.41) and ( 15.791) . and r(-) denotes the PDF for the largest Euclidean 
deviation of the two-dimensional standard Wiener process from the origin over the time interval 
| 0 . 1 ], 


Proof. By using (15.171) and (15.181) . it follows that, for any 1 ^ ^ r and 1 ^ k ^ N, 


where 


n,:= 


■J 3 a-'v,h = 

(5.50) 

[02x2(j-l) I 2 02x(n-2j)] ^ R^^” 

(5.51) 


with Opxg denoting the zero (p x g)-matrix. Since each of the rotation matrices Jj is orthogonal, 
(15.501) implies that 



n j w _fe 


(5.52) 


Note that the vector in (15.201) depends on the auxiliary time variable t G [0,1] in a 

continuous piece-wise linear fashion (it is a linear function of f G [{k — 1)/N, k/N] for any 
1 ^ k ^ N). Therefore, it follows from (15.521) that 


[Y maxi^fc^AT \ Vj6k 

y N 


max 

i6[0,l] 


(5.53) 
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Now, consider the veetor-valued map 


aN := 


/ maxi^kiZN \Vj6k 

y n[ 


l^j^r 


f max imiyjvtl 1 
V6[0-d ’ 


(5.54) 


whose entries are described by (15.531) . For any N eN, this map is representable as the eompo- 
sition Oat = (^ o IFV of the map : Z"- ^ C defined above and a Lipsehitz continuous map 
(y9 : C —)■ on the spaee (15.211) given by 


(p{f) := (max \njf{t)\^ . (5.55) 

/ Isgjsjr 

The preimage of any Jordan measurable set B C belongs to the algebra (15.301) . 

and henee, application of Theorem [16] leads to limAr^.+oo F(a)(f^(i?)) = W{ip~^{B)). The lat¬ 
ter convergence implies (15.491) in view of (15.511) . (15.541) and (15.551) . Here, we have also used 
the property that the projections of the 2r-dimensional standard Wiener proeess onto r pair¬ 
wise orthogonal two-dimensional subspaces are mutually independent two-dimensional stan¬ 
dard Wiener proeesses). ■ 


5.3 An application to the rounded-off planar rotations 


We will now apply the above results to a dynamieal system on the two-dimensional lattiee Z^, 
with which a celebrated problem on the rounded-off planar rotations ||4l[II]l is eoneerned. More 
preeisely, consider the quantized linear (i?*, L)-system, where i?* : —>■ Z^ is the roundoff 

quantizer with 

^*■'( 0 ) = [- 1 / 2 , 1 / 2)2 ( 556 ) 


(see Section 1431) . and 


L := 


cos 9 
sin 6 


— sin 6 
cos 9 


(5.57) 


is the matrix of rotation by angle 9 E (0, 27r). This dynamieal system, whose transition operator 
is given by 

T = R^o L, 


is a partieular ease of the neutral quantized linear systems deseribed in Seetion 15.11 Indeed, 
the matrix L in (15.571) is of the form (I5.1I) - (I5.3I) . with r = 1 and U = R- Note that the 
assumption for the matrix L to have a nondegenerate speetrum is equivalent to 9 ^ n whieh 
holds, for example, if 6^ E (0,7r/2). Furthermore, it is assumed throughout this subseetion that 
the rotation angle 9 belongs to the set 


0 := G (0,7r/2) : the matrix L in fl5.57l) is iteratively nonresonant}; (5.58) 

see Seetion |431 As diseussed in Seetion [5H the set 0 is of full Lebesgue measure, and more¬ 
over, the set (0,7r/2) \ 0 of “pathologieal” values of the rotation angle is eountable. Also, 9/7r 
is irrational for any 9 E Q. 

The following theorem provides corollaries from the corresponding results of the previous 
subsection and Section |4| For its formulation, we note that, in view of (15.561) and (15.571) . the 
mean vector jj, in (14.621) vanishes. 


/i 


udu 


'[- 1 / 2 , 1 / 2)2 


0 , 
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and hence, the supporting system with the transition operator in (14.611) coincides with the origi¬ 
nal linear system. Also note that the matrix <l> in (15.41) and the covariance matrix 4^ in (15.191) for 
the uniform distribution over the square in (15.561) take the form 

$ = ^ = / uu^du = —. (5.59) 

J[-l/2,1/2)2 12 

Theorem 19 Suppose the rotation angle 6 belongs to the set 0 in ( 15.5^1) . Then the quantized 
linear (i?*, L)-system, specified by ( 15.5(51) and ( 15.571) . satisfies the following properties: 

(a) the quantization errors E). given by ( 14.591) are mutually independent and uniformly dis¬ 
tributed over the square [—1/2,1/2)^. In particular, 

( N \ N 

= J]mes25fc 
fc=l / k=l 

holds for any G N and any Jordan measurable subsets Bi,...,B]^ofthe square; 

(b) for any Jordan measurable set B C 


lim F 

A^^ + OO 


X e 



rfix)) G B 


1 



Hi , 

2 au. 


(5.60) 


In particular, for any a ^ 0, 


lim F 

A^^+OO 


X G 



(5.61) 


(c) for any a 0, 


lim F 

A^^ + OO 


X G Z^ 



12 

max 
JV 


jL^x-T'^(x)j > a 


r"-|-oo 


T(u)du, 


where r : M+ —)■ M+ is the PDF of the largest Euclidean deviation of the two-dimensional 
standard Wiener process from the origin over the time interval [0,1]; 

(d) the cardinality of the preimage of a point under the transition operator T takes values 
0,1,2 with the following frequencies 


( fi for A; = 0 

F{xgZ 2 : #T-^(x) = A:} = <^ 1-2/3 for A; = 1 , (5.62) 

y (3 for A; = 2 

where 

/3 := 1 - YiTil?)) = (cos9 + sin 9 - 1)^ (5.63) 

is the frequency of “holes” in the lattice if which are not reachable for the system. 
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Proof. The assertion (a) of the theorem is a corollary from Theorem |9l The assertion (b) is 
established by noting that (15.601) follows from (15.591) and Theorem [iTl The relation (15.611) can 
be obtained from (15.601) by using the fact that the squared Euclidean norm of a two-dimensional 
Gaussian random vector with zero mean and identity covariance matrix has the ^^-distribution 
with two degrees of freedom: 

1 f _a 

— e ^ du = e '2 for all a ^ 0. 

JuGR'^: \u\‘^>a 

The assertion (c) follows from Theorem [TSl The assertion (d) is a corollary from Theorems [13] 
and Us] Indeed, by using (14.801) . it follows that 

F{x e : W-\x) = k} 

= mes2 [ue [-1/2,1/2)^ ; # p|(M + L-^ [-1/2,1/2)2)) ^ 

In combination with (14.911) and (14.921) (see also Fig. [3]), the relation (15.641) leads to (15.621) and 
(15331) . ■ 



Figure 3: Calculation of the frequency F(Z2 \ T{Z'^)) of “holes” in the two-dimensional lattice 
Z2 which are not reachable for the rounded-off planar rotation with a generic angle 9. This 
frequency is equal to the area (cos 9 + sin 6^ — 1)^ of the shaded square. 


Unlike the statements (a)-(c) of Theorem [T9l its assertion (d) involves only the first iterate 
of the transition operator T and remains valid for nonresonant (that is, not necessarily iteratively 
nonresonant) matrices L in (15.571) : see Section |431 Since for the orthogonal (2 x 2)-matrices, 
the nonresonance property is equivalent to irrationality of e*®, then Theorem [T9( dl holds, in 
particular, for 9 = |. In this case, the rational independence of the rows of the corresponding 
matrix 


r 1 0 1 

h 

L 

1 V3 

1-2 2-1 


0 1 

V3 _ 1 
2 2 
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Figure 4: A fragment {—50,..., 50}^ of the lattiee 1?, with the grey pixels depieting the image 
T{I?) of the lattice under the rounded-off planar rotation with angle 0 = ^. The white pixels 
represent the complement consisting of “holes” in the lattice which are not reachable 

for the discretized system. 


can also be verified directly. For this rotation angle, a fragment of the reachability set T{1?) is 
shown in Fig.|4l This is an example of a frequency measurable L“^-quasiperiodic subset of the 
lattice, which is not spatially periodic. According to (15.631) of Theorem [T9{d), the frequency of 
the reachability set is 

F(T(Z2)) = 1 - fcos- + sin- - iV = — = 0.8660... (5.65) 

V 6 6/2 

The relative fraction of reachable points in the moderately large fragment of the lattice under 
consideration is 0.8659..., with the relative error of the theoretical prediction in (15.651) being 
0.015%. 

Also note that, by Theorem [T9(d). the frequency /9 of “holes” for the rounded-off planar 
rotation in (15.631) asymptotically achieves its largest value (\/2 — 1)^ = 0.1716... for generic 
rotation angles 9 approaching For such rotation angles, the transition operator T manifests 
“minimal surjectivity”. 

A comparison of the above discussed theoretical predictions with experimental results can 
also be found in [(H, where the relative fraction of points in sufficiently large rectangular frag¬ 
ments of were calculated (along with their frequencies) for other events relevant to the phase 
portraits of the rounded-off planar rotations. That comparison also demonstrates close proxim¬ 
ity of the numerical experiment and the predictions provided by the frequency-based approach. 
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